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ABSTRACT 


Modern naval warfare has been increasingly dependent upon the 
acoustic silencing of the participants. Constrained viscoelastic layer 
damping of vibrating elements is one method which can be used to meet 
acoustic silencing goals. This paper considers constrained viscoelastic 
layer damping treatments applied to a thick aluminum plate, including 
single layer, double layer, a milled pocket plate, and a milled “floating 
element” configuration. High modal damping values were obtained for 
each damping configuration. The Modal Strain Energy method, using 
finite element analysis to estimate modal loss factors, was investigated for 
use as a tool in constrained viscoelastic layer damping design. A 
comparison of experimentally measured frequency response and modal 
loss factors with those predicted by the modal strain energy method is 
presented to confirm the possible use of the modal strain energy method as 
a design tool. 
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L INTRODUCTION 


A. BACKGROUND 

Modem naval warfare, especially undersea warfare, depends heavily 
on the vessel being acoustically silent. A major source of radiated noise is 
the vibration of shipboard components. The reduction of these vibrations is 
of utmost importance if a ship is to accomplish its mission. One method of 
vibration damping that shows promise in damping over a broad spectrum 
of low frequency vibration is constrained viscoelastic layer damping. The 
constrained viscoelastic layer method uses shear deformation in 
viscoelastic material to absorb and dissipate the vibrational energy of the 
system in question. Unfortunately the design and analysis of such 
constrained viscoelastic layer systems is difficult, due in part to the 
following: 

• The material properties of viscoelastic damping materials vary 
greatly with temperature and frequency. 

• Closed form solutions to the equations of motion exist only for beams 
and plates with simple boundary conditions. 

• The necessity for large amounts of computer storage and CPU time to 
conduct finite element analysis of constrained viscoelastic layer 
damping systems because of a large number of elements and the 
variation of material properties with frequency. 
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A finite element technique developed by Johnson and Kienholz 
known as the Modal Strain Energy (MSE) method uses structural strain 
energy to approximate the damping of a structure with an applied 
constrained viscoelastic damping system [Ref. 1]. Previous work by Maurer 
addressed the effectiveness of the MSE method for two damped plate 
configurations: 1) a simple sandwich configuration, and 2) a plate with a 
milled pocket with damping material inserted and a welded cover plate 
acting as a constraining layer [Ref.2]. Difficulties with the welded cover 
plate were reported because during welding, the cover plate warped and 
delaminated itself from the damping material, resulting in negligible 
damping [Ref. 2]. 

B. PURPOSE 

For certain naval applications the components to be damped will be 
thick in construction and may be exposed to an unfriendly environment. 
This paper addresses the experimental testing and analysis of four thick 
aluminum plates, each with a different constrained viscoelastic layer 
damping treatment. 

This research further investigated the pocket plate configuration and 
also investigated a constrained viscoelastic layer damping system using 
two separate layers of damping material. In addition, a second milled plate 
using a “floating element” in conjunction with constrained layer damping 
was tested. The MSE method is used to analyze each of the structures, and 
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its accuracy and usefulness as a possible design tool are investigated. 
Although the intended use of the MSE method is the design of complex 
damping structures, simple plate geometries were used to facilitate the 
experimental and computational effort. 


3 




II. THEORY 


A. VISCOELASTIC MATERIAL 

Viscoelastic materials of interest for general naval applications are 
polymeric compounds made up of long molecular chains. These molecular 
chains can be strongly, or weakly, linked together, depending on their 
chemical composition and processing. The damping characteristics of 
viscoelastic arise from the deformation and recovery of the polymer 
network. Material properties of a viscoelastic material vary with 
temperature and frequency. As such, the damping characteristics of a 
system will vary as its operating environment changes. [Ref. 3] 

Temperature will have the greatest effect on the material properties 
of damping materials [Ref. 3]. This effect is shown in Figure 2.1, where 
four distinct regions are observed. The lowest temperature region is the 
glassy region where the material’s storage modulus is at its maximum 
value, and the loss factor is at a minimum. In the glassy region the 
modulus decreases slowly with temperature increase, whereas the loss 
factor increases rapidly with temperature. The second region is the 
transition region where the modulus decreases rapidly with increasing 
temperature and the loss factor reaches its peak value. The third region is 
the rubbery region where both the modulus and loss factor are at low values 
and show little variation with temperature. The fourth, and last, region is 
the flow region and characterizes the behavior of some materials, mostly 
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ceramics, at high temperatures. It should be noted that the transition 
region may vary in width from 20 °C up to a width of 200 °C. [Ref. 3] 

The effect of frequency on viscoelastic materials is not as great as that 
of temperature. The modulus of the viscoelastic always increases with 
increasing frequency. The loss factor will initially increase with frequency, 
then peak, and subsequently decrease as frequency increases. A plot of 
storage modulus and loss factor versus frequency is shown in Figure 2.2. It 
should be noted that this plot is over a range of approximately ten decades, 
and hence it becomes obvious that a temperature change of a couple degrees 
will have a much greater effect on damping than a minor change in 
frequency. [Ref. 3] 



Figure 2.1. Variation of viscoelastic material 
properties with temperature [Ref. 3]. 
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Figure 2.2. Variation of viscoelastic material properties 
with frequency and constant temperature [Ref. 3]. 


Linear viscoelastic materials behave in a hysteretic manner under 
cyclic excitation. Therefore, the mechanical behavior of a viscoelastic 
material during steady state vibration is best described by using a complex 
stiffness, k* [Ref. 3]. 

k* = k(l + irj) (2.1) 

where, 

T| = material loss factor 








The use of a complex stiffness then leads to the use of a complex Young’s 
modulus and shear modulus [Ref. 3]. 

E* = E(1 + irt) (2.2) 

G* = Gd + iri) (2.3) 

This concept of the complex modulus is used in subsequent analysis. 

Viscoelastic material properties are commonly displayed using a 
“reduced frequency nomogram.” The reduced frequency nomogram 
displays the variation of the viscoelastic material’s loss factor and modulus 
with temperature and frequency. The “reduced frequency”, f«t. is an 

empirically determined function that accounts for the viscoelastic’s 
temperature and frequency dependence, and allows data for wide range of 
temperature and frequencies to be plotted on the same graph [Ref 4]. The 
reduced frequency nomogram for 3M ISD - 112 is shown in Figure 2.3. To 
find the loss factor and modulus using the nomogram, enter with the 
desired temperature and frequency. Follow the frequency line horizontally 
and the temperature line diagonally down the page until the two intersect. 
Then go vertically up or down to intersect the shear modulus or loss factor 
curves. Finally, read the value of the shear modulus or loss factor 
horizontally from the scale on the left [Ref. 4]. 
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DAMPING PROPERTIES 



8 




























































R CONSTRAINED VISCOELASTIC LAYER DAMPING 

A simple constrained layer damping treatment consists of a base 
layer (the structure to be damped), a damping layer, and the constraining 
layer. This configuration is shown in Figure 2.4 with the thicknesses of the 
damping and constraining layers exaggerated for clarity. 



Figure 2.4. Simple constrained viscoelastic layer configuration. 

The physical mechanism of damping can be explained by referring to 
Figure 2.4. When the base layer is deformed in a mode of vibration, the 
surface away from the neutral axis elongates, stretching the viscoelastic 
material. The top layer, being a stiff elastic material, tends not to elongate, 
and thereby “constrains” the viscoelastic material. Consequently, the cyclic 
motions of vibration induce a cyclic shearing strain in the viscoelastic. 
This cyclic shearing strain, together with its associated hysteresis loop 
cause the vibrational energy to be dissipated as heat. For the constraining 
layer to be effective, its stiffness should not exceed that of the base layer. 
[Ref. 4 & 5] 


9 





C. SYSTEM EQUATIONS OF MOTION 


Continuous systems, such as plates, possess distributed 
characteristics of mass, damping, and stiffness. Classical vibration 
analysis of such systems involves the formation of a mathematical model 
that discretizes the system into a finite number of components in order to 
approximate the total system. Such a formulation results in the following 
equation: 

[M] (x(t)} + [C] {x(t)} + [K] {x(t)} = {F(t)} (2.4) 

where, 

[M] = system mass matrix 
[C] = system damping matrix 
[K] = system stiffness matrix 
(F(t)} = external excitation vector 
(x(t)} = displacement vector 

For an undamped system without excitation, the above equation reduces to 
the eigenvalue problem. 


[M]{x(t)}+[K]{x(t)} = 0 (2.5) 

This equation is then transformed to modal space using the linear 
transformation: 


(x(t)} =[d>] (q(t)} 


(2.6) 






where, 


[4>] = modal matrix 

(q(t)} = modal response vector 

Using this linear transformation the equation of motion can then be solved 
for the undamped modal frequencies and mode shapes. 

To solve for the frequency response of a damped system, the linear 
transformation is applied to equation (2.4): 

[M] M |q(t)) + [C] [<t>] {q(t)J + [K] [<j>] {q(t}} = (F(t)} (2.7) 

Assuming that the damping matrix [C] is proportional to a linear 
combination of the stiffness matrix [K] and mass matrix [M], the damping 
matrix can then be diagonalized using the same linear transformation 
used to diagonalize [K] and [M] in equation (2.5) above. The diagonal terms 
of the damping matrix then become (TijCOj), where Tij equals the modal loss 
factor and coi is the natural frequency of the i^ mode [Ref. 1] . Using this 
approximate diagonal damping matrix results in a system of uncoupled 
modal equations of motion: 

cji(t) + riicq<j(t) + (opqi(t) = fi(t) (2.8) 


where, 

qj(t) = modal acceleration of i^ mode 
q(t) = modal velocity of i^ 1 mode 
qj(t) = modal displacement of i^ mode 
"Hi = modal loss factor of i^ mode 




coj = natural frequency 
fi(t) = modal force in mode 

j = VT 

Assuming that a sinusoidal excitation produces a sinusoidal response, 

(fU)} = (f!ei“* {q(t)} = {Q}ei«* (2.9) 

the response for the i^ mode is then solved to be: 

q . = —J_- (2.10) 

oaf - or + jcDiTtito 

Subsequently, the response of the physical system can be found using: 

{ x (t)} = [<t>]fQ}e><* (2.11) 

D. MODAL STRAIN ENERGY METHOD 

The equations of motion used to define the response of a system with 
viscoelastic materials need a complex eigenvalue analysis. However, the 
actual solution of these equations may be quite difficult. This is especially 
true when the system to be analyzed is comprised of materials whose 
properties vary with both temperature and frequency. Finite element 
techniques are generally used to compute the response of complicated 
systems. However, for the case of varying material properties many time 
consuming and costly runs must be made with the material properties 





changing at each frequency increment [Ref. 1]. In addition to the costly 
analysis of a single design configuration, changes in design options, design 
requirements, of the search for an optimum design can make the expense 
of finite element analysis too great. The development of the Modal Strain 
Energy (MSE) method by Johnson and Kienholz, however, makes the finite 
element analysis of complex viscoelastically damped structures a viable 
option [Ref. 1], 

The MSE method assumes that a damped structure can be 
represented by the normal modes of the associated undamped system if 
appropriate damping terms are inserted into the uncoupled modal 
equations of motion (eqn 2.8) [Ref. 1]. The MSE method further assumes 
that the damping matrix of equation (2.4) can be diagonalized by the same 
real modal matrix that diagonalizes the system mass and stiffness 
matrices of equation (2.4). 

The modal loss factors of equation (2.8) are calculated using the 
undamped mode shapes and material loss factors for each material. Since 
the material loss factors of the structure to be damped and the constraining 
layers are quite small compared to those of common viscoelastic cores, the 
modal loss factors can be estimated using the following equation [Ref. 1]: 

TJ (r) = ~TT (2.12) 


where, 

Tl (r) = system loss factor in the r^ mode 





Ti ( v r) = material loss factor of viscoelastic core at the r^ 
natural frequency 

Vy r) = strain energy in the viscoelastic core at r^ mode 
V (r) = strain energy in the entire structure at r^ mode 

The modal strain energies can be obtained for finite element analysis, and 
are a standard output option of the NASTRAN finite element code [Ref. 7], 
The modal frequency response of the structure is then calculated 
using the modal loss factors found in equation (2.12). When computing the 
modal frequency response of a damped structure, the modal properties in 
the system matrices are assumed to be constant. However, viscoelastic 
materials have storage moduli which are frequency dependent. To account 
for this frequency dependence, Johnson and Kienholz devised the following 
correction factor to be applied to the modal loss factors calculated in 
equation (2.12) [Ref. 8]. 



(2.13) 


where, 

= corrected modal loss factor at the r^ 1 mode 
G 2 (f r ) = viscoelastic shear modulus at the r** 1 modal 
frequency 

G 2 ,ref = reference viscoelastic shear modulus used in the 
frequency response calculation 




IIL DESIGN OF DAMPED PLATES 


A. GENERAL SPECIMEN CONFIGURATIONS 

For experimental testing and analysis purposes, four different, yet 
related, constrained layer damping treatments were selected in addition to 
an undamped “reference” plate. Two of the damping configurations were 
simple sandwich treatments consisting of one and two viscoelastic layers 
respectively. Another treatment was the “pocket plate” which was 
previously investigated by Maurer [Ref. 2]. The pocket plate was made from 
a solid plate which was then milled to accept damping material and a cover 
plate. The final damping treatment was a “floating element” configuration 
consisting of a solid plate milled to accept a double layer of damping 
material and a welded cover plate. Section views of these damping 
configurations are shown in Figure 3.1. 

The purpose of the pocket plate is to protect the viscoelastic material 
from materials such as oil and salt water, which may harm the 
viscoelastic. Since previous attempts at using a welded cover plate were 
unsuccessful due to the heat of welding causing a delamination between 
the damping material and cover plate [Ref. 2], it was decided to move the 
viscoelastic material from the welding point by recessing it into a shallow 
pocket of its own as shown in Figure 3.1. 

The “floating element” concept evolved because the welded cover plate 
of the pocket plate configuration does not produce the damping reaction that 
a true constraining layer would provide. If the cover plate is welded to the 




surrounding structure it cannot deform in bending as much as an 
unwelded constraining layer, thereby causing a reduction in the damping 
capability of the system. By using a piece of metal in the milled pocket with 
dimensions slightly smaller than the surrounding pocket, along with two 
layers of viscoelastic and a welded cover plate, a true constraining layer 
effect should be obtained. This configuration is shown in Figure 3.1. 




Single layer configuration 


Double layer configuration 


I rr 7 > > ? s 7 / / ; j mr 


Milled “pocket plate” 




Plate with “floating element” 


Figure 3.1. Generalized damping treatment configurations. 


In order to approximate a possible system to be damped, a plate with 
large dimensions was selected. The dimensions of the plates used for the 
damping treatments are 114.3 cm (45 in) in length and 38.1 cm (15 in) in 
width. In addition, possible naval applications for this type of damping 
would probably consist of thick plate members. For this reason it was also 
decided that thick plates would be used for the damping treatments. All the 







base plates and constraining layers were made of a standard 6061-T6 
aluminum alloy. The design and selection of viscoelastic layer, base layer, 
and constraining layer thicknesses is discussed in the following sections. 

B. DESIGN OF THE SINGLE DAMPING LAYER CONFIGURATION 

In an attempt to approximate system loss factors and hence 
determine viscoelastic and constraining layer thicknesses for maximum 
damping, a method de/eloped by Nashif [Ref. 9] based on an analysis of 
simple sandwich damping systems by Ross, Kerwin, and Ungar [Ref. 10] 
was used. The Ross - Kerwin - Ungar (RKU) equations are base of the 
analysis of the simple sandwich system shown in Figure 3.2. 


CONSTRAINING LAYER 


DAMPING LAYER 


BASE LAYER 



Figure 3.2. Elements of a simple sandwich damping system. 





To find the loss factors of the damped system, the flexural rigidity of 
the system must first be determined. For the above system, the flexural 
rigidity, El, is written as [Ref. 9]: 


El = Eill + + EsMl + E Hl D 2 

12 12 12 

+ E 2 H 2 (H 21 - D ) 2 + E 3 Ha(H3i - D ) 2 


E^l + E^Ha, Dl + EaHsfHa, D)]Sai_J> 

L 12 2 J 1 + g 

(3.1) 

where, 


E 2 H 2 (h 2 i - MSL) + g(E 2 H 2 H 21 + E 3 H 3 H 31 ) 

D = ' 2 ' 

EiHj + + gfEjHx + E 2 H 2 + E 3 H 3 ) 

(3.2) 

H 31 = Hl ' + Ha t H 2 

(3.3) 

h 21 = h . + 2 H 2 

(3.4) 

g= ° 2 2 

e 3 h 3 h 2 k 2 

(3.5) 


E = Young’s modulus 
G = Shear modulus 
I = second moment of area 
H = thickness of member 
K 2 = wave number 

Subscripts refer to the layers labeled in Figure 3.1 
No subscript refers to the composite system 




For a simply supported plate the wave numbers and modal frequencies are 
found using [Ref.9]: 


— Knm 


EH 3 g c 
12(l - v 2 )Hp 


(3.6) 



where, 

a = semi—wave length of the plate 
b = semi-wave width of the plate 
v = Poisson’s ratio of the composite plate 
p = density of the composite plate 
g c = gravitational constant 


(3.7) 


To introduce damping into the equations it is necessary to use the 
complex modulus concept expressed in Section II. Substituting the 
appropriate expressions for the complex shear and Young’s modulus into 
equations (3.1), and assuming that damping in the base layer, (ri 1 ), is 

small, and that the extensional stiffness of the damping layer is small 
(since E 2 « E^ and E 2 « E 3 ), the following expressions can be arrived at 

[Ref.9]: 


EH 3 = ExH? + E 3 Hj] + —12— (a - p - 8)r E 

c 2 + d 2 


1 9 


(3.8) 





(3.9) 


EH 3 !! = E 3 H 3 TI 3 + — 12 —(a -P-3U 

c 2 + 


where, 

a = gEiHxEaHsHiilcfl -TfeTts) + $t \2 + t\ 3 ) + j[c(rj 2 + tls)]} (3.10) 
p = E 1 H 1 E 2 H 2 H 31 [c + dri 2 + j(cT ]2 - d)] (3.11) 


d = 2gE 2 H2E3H 3 H2iH3i 


c(l -2Tl 2 Tl3-TlD + d( 2 Tl 2 + Tl3-Tl2n3) 1 
{c(2r\2 + Tl3 - TI 2 TI 3 ) - d(l - 2T12T13 - Til)]) 


(3.12) 


c = EiHi(l + g) + gE 3 H 3 (l - TI 2 TI 3 ) (3.13) 

d = gExHiTi2 + gE 3 H3(Ti2 + T1 3 ) (3.14) 

j = (3.15) 


These equations were then applied to estimate the loss factors of the simple 
three-layer sandwich plate. The equations can be simplified by assuming 
that damping in the constraining layer ( 1 ) 3 ) is negligible [Ref.9]. 

Since the boundary conditions for the plate used in this research 
(free-free-free-free) do not correspond to the simply supported conditions on 
which equation (3.7) is based, modal frequencies for the free-free case were 
estimated using results from finite element analysis of a free-free plate and 
equation (3.6). The natural frequencies of an undamped plate with the 
dimensions previously given, and a thickness of 1.91 cm (0.75 in) were 
found using a normal mode extraction in NASTRAN. By substituting these 
modal frequencies into equation (3.6) and estimate of the wave parameter 
for each mode, K^m, was obtained. Then, by using an iterative procedure 
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outlined in Reference [9], modal loss factors were estimated for different 
layer thicknesses over a temperature range of 0.0 °C to 37.8 °C (30 °F to 
100 °F). 

The previous equations are easily programmed to compute loss 
factors for a wide variety of conditions. The variation of viscoelastic 
material properties with temperature and frequency was accounted for 
using a curve-fit to the reduced frequency nomogram developed by Drake 
[Ref. 11]. The material data for the following curve-fit equations is from the 
University of Dayton Research Institute [Ref. 12]. 


2 


logi 0 {M) = logio(ML) + 



(3.16) 


log 10 (ETA) = log 10 (ETAFROL) 

+ lc[A(SL + SH) + (SL - SH/l - Vl + A 2 )] ( ' } 

jL 

logiotFR) ■ logio(F) - (6 25 ^ j ^T 0 ( <318> 


A logi o(FR) - logt o(FROL) 

A " C 

where, 

M = viscoelastic modulus 

ETA = viscoelastic material loss factor 

FR = reduced frequency (Hz) 

F = frequency (Hz) 

T = temperature (°F) 


(3.19) 





and, 


T0 = 40 °C (104 °F) 

FROM = 2.0x10 4 Hz 

MROM = 4.75x10 6 Pa (688.94 psi) 

n = 0.275 

ML = 6.0xl0 4 Pa (8.7 psi) 

ETAFROL = 1.08 
SL = 0.45 
SH = -0.55 
FROL = 5000 Hz 
C = 2.5 

In addition to the above constants, a Poisson’s ratio of 0.49 and a density of 
0.909 gram per cubic centimeter (0.035 Ibm/in®) was used for ISD-112 
[Ref.12], The following material properties were used for 6061-T6 
aluminum [Ref.13]. 

E = 70 GPa (10x10® psi) 
v = 0.33 

p = 2.7 gm/cm® (0.0968 lbm/in®) 

Using the previous equations and material properties, a computer 
program was written to compute estimated modal loss factors and 
frequencies for a variety of base layer, viscoelastic layer, and constraining 
layer thicknesses. A listing of this program appears in Appendix A. 
Modal loss factors were computed for base layer thicknesses of 9.53 mm 
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(0.375 in) to 19.05 mm (0.75 in) in 3.18 mm (0.125 in) increments. For each 
base layer thickness, the viscoelastic thickness was varied from 0.38 mm 
(0.015 in) to 1.52 mm (0.060 in) in 0.38 mm increments, and the 
constraining layer thickness was varied from 1.59 mm (0.0625 in) to 
6.35 mm (0.25 in) in 1.59 mm increments. In addition, loss factors were 
also computed for a viscoelastic thickness of 0.127 mm (0.005 in). From the 
results of the analysis, a carpet plot [Ref. 3], was made for each of the base 
layer conditions. The carpet plot reflects, for the the first mode, the 
maximum loss factor and its corresponding temperature for each 
viscoelastic layer/constraining layer thickness configuration. The carpet 
plot for a base layer thickness of 12.7 mm (0.50 in) is shown in Figure 3.3. 
Based on the carpet plots and a desire for maximum damping, as well as a 
system which could be moved easily, a base layer thickness of 12.7 mm 
(0.50 in), a viscoelastic thickness of 0.38 mm (0.015 in), and a constraining 
layer thickness of 6.35 mm (0.25 in) was selected. The total system 
thickness was approximately 19.05 mm (0.75 in). This total system 
thickness would be maintained for all subsequent damping configurations. 

To maintain continuity between damping systems, the milled “pocket 
plate” was given a viscoelastic thickness of 0.38 mm and a cover plate 
thickness of 6.35 mm for a total system thickness of 19.05 mm (0.75 in). In 
order to keep the heat of welding away from the viscoelastic material, the 
ISD-112 was recessed into a shallow pocket as indicated previously in 
Figure 3.1, and as shown in the pocket plate system arrangement in 
Figure 3.4. Detail drawings of the pocket plate are shown in Appendix C. 
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Maximum System Loss Factor 



Figure 3.3. Carpet plot for a base layer thickness, HI = 12.7 mm. 
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C. DESIGN OF THE DOUBLE DAMPING LAYER 


The design of the double layer damping system was accomplished 
using a modification of the RKU analysis used in the previous section. The 
RKU equations are used by working from the top layer of the damping 
system down towards the base layer. As shown in Figure 3.5, the H3' 
constraining layer along with the H2' viscoelastic layer are combined with 
the HI' layer to form a three-layer system. Using the RKU equations, the 
stiffness of this system is computed and considered to be the equivalent 
stiffness of the top three layers of the total constrained layer damping 
system. The top three layers were then considered as a single layer with 
the equivalent stiffness previously calculated, and the RKU equations were 
again applied to compute estimated modal loss factors for the entire double 
layer damping system. [Ref.3] 


CONSTRAINING LAYER 
VISCOELASTIC 


CONSTRAINING LAYER 
VISCOELASTIC 


BASE STRUCTURE 


H3' 

H2’ 

HI' 

H2 


HI 


Figure 3.5. General configuration of the double 
constrained layer system. 
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To maintain continuity among all the damping configurations, a 
base layer thickness of 12.77 mm (0.50 in) was chosen, and a total system 
thickness of 19.05 mm (0.75 in) was maintained. A design for high 
damping was then selected by computing modal loss factors for the 
constraining and damping layer thickness combinations shown in Table 
3.1. The estimated modal loss factors for the first mode of vibration in each 
configuration are plotted as shown in Figure 3.6. A listing of the program 
used to compute the loss factors is in Appendix B. 

From the data presented in Figure 3.6, viscoelastic thickness of 
0.38 mm (0.015 in) and constraining layer thicknesses of 3.18 mm (0.125 in) 
were selected for the double layer configuration. This particular 
configuration estimates high damping over a wider temperature range 
than the other thickness combinations. 


TABLE 3 X THICKNESSES USED IN CALCULATION OF DOUBLE 
LAYER MODAL LOSS FACTORS. 


Configuration 

1 

2 

2 

4 

5 

6 

H2 (mm) 

0.38 

0.38 

0.38 

0.76 

0.38 

1.14 

HI' (mm) 

3.18 

2.38 

3.18 

2.38 

2.38 

2.38 

H2' (mm) 

0.38 

0.38 

0.38 

0.76 

1.14 

0.38 

H3' (mm) 

3.18 

3.18 

2.38 

2.38 

2.38 

2.38 
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Figure 3.6. Modal loss factors for double constrained layer 
damping configuration. 

The milled “floating element” plate uses the same viscoelastic and 
constraining layer thicknesses as the simple two-layer configuration. In a 
design similar to that of the pocket plate, the floating element and both 
layers of viscoelastic are recessed into a milled opening as shown 
previously in Figure 3.1 and further described in the floating element 
system configuration of Figure 3.7. 
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Figure 3.7. Floating element system configuration. 
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IV. EXPERIMENTAL RESULTS 


A. TESTING ARRANGEMENT 

Experimental testing was performed on each of the four damping 
configurations and the undamped reference plate. In order to approximate 
the free-free-free-free boundary condition, each plate was suspended from 
the roof of the testing chamber using elastic cords as shown in Figure 4.1. 
All of the tests were performed in a temperature controlled environmental 
chamber which enabled temperatures to be maintained within ±1 °C. The 
primary component and user interface was the Hewlett-Packard (HP) 
3562A Dynamic Signal Analyzer (DSA). The HP-3562A was used to provide 
a swept sine signal to a vibration generator, and analyzed the returning 
data signals. The HP-3562A was used to compute the frequency response 
and coherence over a range of 50 Hz to 1050 Hz using the discrete Fourier 
Transform in the swept sine mode. Ten averages were performed at each 
data point using a frequency resolution of 625 mHz per step. The source 
level output to the vibration generator was set at 1.5 volts. 

A schematic of the testing apparatus is shown in Figure 4.2. Swept 
sine source signals were fed from the output jack of the HP-3562A to a 
Wilcoxon F3 vibration generator via the piezoelectric output of a Wilcoxon 
PA7C power amplifier. The vibration generator was mounted 73.48 cm 
(28.93 in) from one end, and 12.7 cm (5.0 in) from the front edge of each 
specimen as shown in Figure 4.3. An integral force transducer was 
mounted in the base of the vibration generator to measure the force input to 
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the plate. This force signal was then fed to input channel one of the DSA 
via a PCB 462-A charge amplifier. Plate accelerations were recorded at 
various points using a PCB 303A-03 accelerometer as shown in Figure 4.3. 
Acceleration data was fed to input channel two of the DSA via a PCB 482A05 
power supply. Frequency response and coherence data was then recorded 
on disk for further analysis. 

Temperatures within the testing chamber were maintained using a 
NESLAB RTE-8 refrigerated circulating bath which pumped fluid through 
a small heat exchanger in the testing chamber as shown in Figure 4.2. In 
order to accurately monitor the temperature of the plates, a small 
thermocouple was inserted in the base of each plate. 


SUPPORTING BEAM 



Figure 4.1. Configuration of plate in testing chamber. 







Figure 4.2. Experimental testing system schematic. 
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Figure 4.3. View of shaker and accelerometer locations. 


33 








B. TESTING PROCEDURE 


1. Undamped reference plate 

An undamped, reference, frequency response measurement 
was made at a temperature of 15.6 °C (60 °F) to set a standard response by 
which to measure the effectiveness of the damping treatments. The 
undamped frequency response was recorded over a frequency range of 
50 -1050 Hz using a resolution of 625 mHz per point in the DSA. 

2. Damped plate measurements 

Frequency response measurements of the damped plates were 
made at a temperature of 15.6 °C (60 °F) at several nodes on the plates in 
order to capture the damped response of as many modes as possible. A 
representation of these nodes is shown in Figure 4.3. Responses were 
recorded over a frequency range of 50 -1050 Hz with a resolution of 
625 mHz per point in the DSA. Zoom measurements were also made to 
capture better data for certain modes. Modal loss factors were then 
estimated from the frequency response and coherence measurements 
using a curve-fitting technique described in Reference [14]. 

C. SINGLE DAMPING LAYER RESULTS 

The single damping layer treatment was tested at 4.44 °C (40 °F), 

15.6 °C (60 °F), and 26.7 °C (80 °F) so that the effects of temperature on the 
damping treatment could be determined. A plot of the single damping 
layer frequency response at 15.6 °C is shown in Figure 4.4. The single layer 
damping treatment resulted in high damping with modal loss factors 
ranging from 0.223 at 53.4 Hz to 0.091 at 876.6 Hz. Due to the coupling of 
modes, loss factors for all modes were not measured. The frequency 
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response of the damped plate is characterized by a frequency shift to the left 
and a smoothing of the frequency response when compared to the 
undamped reference plate. The single layer treatment was especially 
effective at reducing the frequency response of a mode cluster between 650 
and 950 Hz. The frequency band of this cluster was shifted approximately 
200 Hz with the amplitudes of the responses of the modes being 
dramatically reduced. The single layer treatment was also effective at 
reducing the amplitude of the response peaks over the entire spectrum of 
measurement. On average, the highest peaks of the frequency response in 
the undamped condition were reduced by 25 decibels, a reduction of 17.8 
times. 

The effect of temperature on the damping was quite pronounced as 
shown in Figure 4.5. As the testing temperature was decreased, the 
viscoelastic layer became stiffer and damping levels were increased. A 
comparative listing of the loss factors at different temperature is in Table 
4.1 and a plot of the modal loss factors is shown in Figure 4.6. Figure 4.5 
shows the trend of increased damping with temperature decrease, and a 
corresponding shift of modal frequencies to the right as the viscoelastic 
becomes stiffer. These changes are especially discernible at the lower 
frequencies. 







Fxd y 50 Hz SINGLE 15.6 C 1.05k 


Figure 4.4. Frequency response of the single damping layer at 16.6 °C. 
t-: undamped response; -: damped response] 



Figure 4.5. Frequency response of the single damping layer configuration 

at different temperatures. 

[-: 4.44 °C,-: 15.6 °C,-: 26.7 °C] 
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TABLE 4.L MODAL LOSS FACTORS FOR THE SINGLE DAMPING 
LAYER AT DIFFERENT TEMPERATURE. 


4.44 °C (40 °F) 15.6 °C (60 °F) 26.7 °C (80 °F) 


££Hzl 

a 

UEzl 

a 

nm 

a 

64.3 

0.217 

53.7 

0.223 

49.0 

0117 

102.8 

0.223 

89.1 

0145 

83.5 

0.089 

157.4 

0.203 

138.9 

0172 

129.6 

0.082 

208.4 

0.183 

188.1 

0184 

177.9 

0.072 

327.5 

0.158 

301.6 

0120 

243.9 

0.062 

440.1 

0.144 

424.8 

0109 

409.4 

0.053 

474.9 

0.205 

476.8 

0111 

450.3 

0.066 

653.1 

0.079 

608.1 

0.067 

552.9 

0.068 

768.0 

0135 

644.2 

0.069 

586.7 

0.052 

852.0 

0.093 

722.5 

0130 

629.3 

0.042 

919.0 

0104 

817.6 

0.081 

681.4 

0.096 

953.5 

0.073 

876.6 

0.091 

790.1 

0.049 


858.2 0.056 

978.1 0.069 



0 200 400 600 800 1000 


Frequency (Hz) 

Figure 4.6. Comparison of modal loss factors for the single layer 
configuration at different temperatures. 
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D. DOUBLE DAMPING LAYER RESULTS 


The double layer damping configuration was also tested at 4.44 °C, 
15.6 °C, and 26.7 °C. The frequency response of this configuration at 15.6 °C 
as compared to the undamped reference plate is shown in Figure 4.7. 
Damping in the double layer configuration is also high, with modal loss 
factors ranging form 0.301 at 53.3 Hz to 0.107 at 832.4 Hz. Due to modal 
coupling loss factors for all modes were not measured. The frequency 
response of the two-layer configuration is also characterized by a dramatic 
reduction in response amplitude and a frequency shift to the left. The peak 
undamped responses were reduced by an average of 27 decibels, or a 
reduction of 22.4 times from the reference condition. 

The effect of temperature on the double layer damping treatment is 
shown in Figure 4.8. As with the single layer case, damping in the double 
layer configuration increased with a decrease in temperature. This 
configuration also shows the shift of modal frequencies to the right as 
temperature decreases and the viscoelastic becomes stiffer. Modal loss 
factors for the double layer configuration are listed in Table 4.2, and are 
plotted for comparison in Figure 4.9. 

To compare the effectiveness of the single layer and double layer 
configurations their frequency responses at 15.6 °C are plotted in Figure 
4.10 with a plot comparing modal loss factors in Figure 4.11. The two 
responses are quite similar although the double layer configuration does 
show an increase of approximately 22 percent in modal loss factor. 
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FxdXY 50 Hz DOUBLE LAYER 15.6 C 1.05k 


Figure 4.7. Frequency response of the double layer at 15.6 °C. 
[-: undamped response ,-: damped response] 



FxdXY 50 Hz A . A A / 15.6 / 26.7 C 1.05k 

Figure 4.8. Frequency response comparison of the double layer 


configuration at different temperatures. 

[-: 4.44 °C,-: 15.6 °C ,-: 26.7 °C] 
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TABLE 42. MODAL LOSS FACTORS FOR THE DOUBLE LAYER 
CONFIGURATION AT DIFFERENT TEMPERATURE 


4.44 °C (40 °F) 


nm 

n 

51.9 

0.273 

122.5 

0.224 

212.7 

0.187 

290.5 

0.197 

382.7 

0.174 

494.2 

0.144 

558.2 

0.169 

602.2 

0.198 

642.4 

0.157 

752.8 

0168 

815.4 

0159 

896.3 

0156 

968.3 

0130 


15.6 °C (60 °F) 


UHzl 

21 

55.3 

0.301 

871 

0.215 

142.1 

0.217 

190.5 

0.212 

297.8 

0.139 

366.8 

0.154 

419.3 

0.125 

441.2 

0100 

618.7 

0.098 

680.9 

0.096 

715.5 

0.060 

832.4 

0107 


26.7 °C (80 °F) 


f (Hz) 

n 

49.7 

0202 

81.9 

0.188 

127.6 

0156 

172.4 

0.117 

278.2 

0.097 

317.7 

0.077 

396.1 

0.070 

428.5 

0.077 

602.0 

0.067 

638.4 

0.072 

738.0 

0.071 

802.6 

0.077 

8461 

0.050 



Figure 4.9. Comparison of modal loss factors for the double layer at 

different temperature. 
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E. POCKET PLATE RESULTS 


The milled pocket plate was constructed as previously shown in 
Figure 3.3. The constraining layer, or cover plate, was welded in place 
using tack welds in an attempt to keep the damping material away from the 
heat of welding, and the cover plate from warping, instead of using a 
continuous weld bead as was previously attempted [Ref. 2]. The cover plate 
was welded to the base at the comers, at the midpoint of the short side and 
at three equally spaced locations along the long dimension as shown in 
Figure 4.12. Following welding the plate was tested to ensure that the 
viscoelastic had not been damaged by the heat of welding. 



Figure 4.12. Location of tack welds on the cover plate 
of the pocket plate configuration. 

The pocket plate was tested at 15.6 °C (60 °F) and the frequency 
response is shown in Figure 4.13. The response indicates that the 
viscoelastic layer was not damaged by welding and that good damping was 
attained. Modal loss factors ranged firom 0.067 at 62.1 Hz to 0.090 at 923 Hz. 
Although damping is good, it is approximately half that of the single layer 
configuration. One reason for this is that the viscoelastic material does not 




completely cover the base structure. Another reason is the presence of the 
welded cover plate. Due to the welded conditions the cover plate cannot 
induce shear deformation in the viscoelastic layer as well as a true 
constraining layer, and therefore produces less damping than the single 
layer configuration. The effects of the welded cover plate are especially felt 
in modes below 300 Hz where the frequency response is quite peaked. The 
response curve becomes more rounded and the effects of the damping layer 
are seen as frequency increases. 

The modal loss factors for the pocket plate are listed in Table 4.3 and 
are plotted in Figure 4.14. Even though the damping is less than the single 
layer, the plate is still adequately damped as shown in the frequency 
response plot in Figure 4.13. In this configuration the modal loss factors 
remained relatively constant throughout the testing spectrum. The 
increase in modal loss factor values above 800 Hz is due primarily to modal 
coupling, and the measured modal loss factors in this range are not 
reliable. 








Fxd V 50 Hz POCKET PLATE 15.6 C 1.05k 


Figure 4.13. Frequency response of the pocket plate at 15.6 °C (60 °F). 
[-: reference plate ,-: pocket plate ] 
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Modal Loss Factors 


TABLE 43 MODAL LOSS FACTORS FOR THE 
POCKET PLATE CONFIGURATION. 


mizi 

21 

62.1 

0.067 

93.4 

0.042 

142.0 

0.056 

194.9 

0.060 

308.9 

0.051 

437.5 

0.041 

495.6 

0.036 


HHz) 

21 

565.3 

0.043 

625.7 

0.047 

643.2 

0.049 

687.0 

0.044 

841.9 

0.056 

891.5 

0.077 

923.0 

0.090 



Figure 4.14. Modal loss factors versus frequency for 
the pocket plate configuration. 




F. 


FLOATING ELEMENT RESULTS 


The milled “floating element” configuration was constructed as 
previously shown in Figure 3.7. The cover plate was welded in a fashion 
similar to the pocket plate as shown in Figure 4.12. The center 
constraining layer, or floating element, was made slightly smaller than the 
surrounding structure thus allowing the floating element to act as a “true” 
constraining layer. 

The frequency response of the floating element configuration at 
15.6 °C (60 °F) is shown in Figure 4.15. The floating element is quite 
effective as the response shows a good reduction in peak modal response. 
Measured modal loss factors range from 0.089 at 66 Hz to 0.064 at 935 Hz. A 
listing of measured modal loss facto rs is in Table 4.4 and are plotted in 
Figure 4.16. As with the previous cases, the frequency response of the 
floating element configuration is characterized by a frequency shift to the 
left and a smoothing of the response as frequency increases. 

In a comparison of the pocket plate and floating element 
configurations, the frequency responses are plotted in Figure 4.17. A 
comparison of modal loss factors for the ‘wo configurations is shown in 
Figure 4.18. The two frequency response plots are similar, however, the 
frequency response of the floating element configuration is more rounded 
than that of the pocket plate. The major difference between the two 
configurations is seen in Figure 4.18. Modal loss factors for the floating 
element show an average increase of 25 percent over those of the pocket 
plate. Reasons for this increase are the added constraining effect of the 
floating element and additional layer of damping material present 
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Modal Loss Factors 


TABLE 4.4 MEASURED MODAL LOSS FACTORS FOR THE FLOATING 

ELEMENT CONFIGURATION. 


mi 

n 

f (Hz) 

n 

66.0 

0.089 

508.1 

0.089 

104.9 

0.040 

540.6 

0.115 

154.0 

0.094 

618.7 

0.121 

211.4 

0.058 

656.7 

0.100 

272.8 

0.063 

708.6 

0.088 

323.5 

0.075 

857.6 

0.089 

435.8 

0.090 

935.4 

0.064 


0.14 
0.12 
0.10 
0.03 
0.06 
0.04 
0.02 

0 200 400 600 800 1000 

Frequency (Hz) 

Figure 4.16. Modal loss factors for the floating element configuration. 
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FREQ RESP 



Figure 4.17. Frequency response comparison for the pocket plate and 
floating element configurations. 

[-: pocket plate ,-: floating element ] 



Frequency (Hz) 

Figure 4.18. Modal loss factor comparison for the pocket plate and floating 

element configurations. 
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V. FINITE ELEMENT RESULTS 


A. UNDAMPED REFERENCE PLATE 

The first step in the finite element analysis procedure was to model 
and analyze the undamped reference plate for its modal frequencies and 
frequency response. The finite element model was generated using 
PATRAN, a computer aided interactive graphics program developed by 
PDA Engineering. PATRAN is widely used for conceptual, preliminary, 
and detailed design and analysis of complex systems. One of PATRAN’s 
major advantages is in the interactive construction of finite element models 
for use by MSC/NASTRAN, and its ability to display MSC/NASTRAN 
results in an easily understood graphic format [Ref. 15]. 

The reference plate was modeled using 84 plate (QUAD4) elements as 
shown in Figure 5.1. The QUAD4 element is an isoparametric element 
with four nodes, one at each comer of the element [Ref. 7]. A normal mode 
extraction was then performed in order to compare numerical results with 
experimentally determined modal frequencies. Once satisfactory 
agreement between the modal frequencies calculated in NASTRAN and 
those obtained experimentally was attained, a direct frequency response 
calculation was performed in NASTRAN. This frequency response was 
then used as the reference response for further finite element models 
incorporating viscoelastic damping treatments. The excitation for the 
frequency response calculation was a sinusoidal force with an amplitude of 
1.0 applied at node 66. The response point was node 58 as shown in Figure 
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5.1. These two nodes correspond to the points on the plates used in the 
experimental portion of the research where the vibration generator and 
accelerometer were attached. 



Figure 5.1. Finite element model of the undamped reference plate. 

B. SINGLE DAMPING LAYER 

The modeling of the single constrained layer damping system was 
done using techniques described by Job "on and Kienholz [Ref. 1]. As 
shown in Figure 5.2, the viscoelastic layer is modeled using solid (HEXA) 
elements, while the base layer and constraining layer were modeled using 
QUAD4 elements. The HEXA element is a solid, isoparametric element 
having eight nodes, one at each corner of the element with three 
translational degrees of freedom at each node [Ref. 7]. The use of solid 
elements for the viscoelastic layer allows the strain energy due to shearing 
to be adequately represented. Plate elements are used in the base layer and 
constraining layer because of their ability to account for stretching and 










bending deformations. The plate element allows its nodes to be offset from 
the plate’s center to the surface of the plate, coincident with the corner 
nodes of the solid viscoelastic elements [Ref. 1]. Thus, the single 
constrained layer system was modeled using only two layers of nodes, a 
simple process in PATRAN. For the single layer configuration, a model 
having 84 elements per layer, and an element meshing scheme shown in 
Figure 5.2 was used. 



Figure 5.2. Finite element representation of the 
single constrained layer. 





Once the damped plate had been modeled, normal mode extractions 
were made using MSC/NASTRAN. Five separate runs were conducted 
using the material properties of ISD-112 at 50, 200, 500, 800, and 1000 Hz 
and at a temperature of 15.6 °C (60 °F). In addition to the modal 
frequencies, the strain energy in the viscoelastic elements and the entire 
model were output from NASTRAN. 

Since the shear modulus of a viscoelastic material changes with 
frequency, it was necessary to estimate the actual modal frequencies of the 
damped plate using an interpolation procedure outlined by Johnson and 
Kienholz [Ref. 16]. The first step of the interpolation process was to plot the 
shear modulus of ISD-112 versus frequency from 5 to 1000 Hz. Then, for 
the first mode, using NASTRAN results based on ISD-112 material 
properties a 50 Hz, the first modal frequency predicted by NASTRAN and 
the corresponding shear modulus were plotted. The same was then done 
using the first natural frequency predicted by normal mode extraction 
based on ISD-112 material properties at frequencies of 500 and 1000 Hz. A 
curve was then passed through these three points. The point where the 
NASTRAN modal frequencies for the first mode intersected the ISD-112 
shear modulus curve was taken to be the interpolated modal frequency of 
the single damping layer configuration. A plot of the intersection of these 
two curves for the first and second modes is shown in Figure 5.3. This 
interpolation process was then repeated for each mode through 1000 Hz. 





Figure 5.3. Interpolation of the first and second modal frequencies 
for the single constraining layer configuration. 

Once the interpolated modal frequencies were found, the modal 
strain energy equations (2.14) and (2.15) were used to compute modal loss 
factors for the single layer configuration. A set of modal loss factors was 
computed based on the modal strain energies computed using viscoelastic 
properties at reference frequencies of 50 ,200, 500, and 800 Hz. A set of 
composite modal loss factors for the modal frequencies near these reference 
frequencies was then selected. The resulting modal loss factors are shown 
in Table 5.1 and are plotted versus frequency in Figure 5.4. As seen in 
Figure 5.4, the modal strain energy method is predicting high damping for 
this configuration with an average modal loss factor of 0.195. 
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Modal Loss Factor 


TABLE 5.L MODAL LOSS FACTORS FOR THE SINGLE 
LAYER CONFIGURATION AS COMPUTED USING 
THE MODAL STRAIN ENERGY METHOD 
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0.180 

322 

0.214 

827 

0.134 

447 

0.200 

869 

0.152 

463 

0.108 

883 

0.115 

483 

0.258 





Figure 5.4. Estimated modal loss factors for the single layer configuration 
with the curve fit used for the MSC/NASTRAN damping table. 






Using the set of composite modal loss factors, the modal frequency 
response of the damped plate was computed using MSC/NASTRAN. Modal 
damping was introduced to the model using the SDAMP option in the Case 
Control Deck and the TABDMP1 damping table in the Bulk Data Deck as 
described in the MSC/NASTRAN Handbook for Dynamic Analysis [Ref. 17]. 
Since NASTRAN uses a linear interpolation between points in the damping 
table to describe the modal damping in the model [Ref. 7], a simple curve fit 
was applied to the set of composite modal loss factors as shown in 
Figure 5.4. Points from this curve fit were then used in the NASTRAN 
damping table. To compute the modal frequency response, a unit excitation 
force was applied at the same node as the undamped plate, and the node 
used for the response was also the same as the undamped plate. 

The results of the modal frequency response calculations are shown 
in Figure 5.5. The dashed line represents the undamped reference plate, 
and the solid line represents the modal frequency response of the single 
layer configuration. Material properties at 200 Hz were used for the 
ISD-112 damping material. The first thirty modes were used in the modal 
summation for the response. A listing of the MSC/NASTRAN data deck 
used to compute the modal frequency response is in Appendix D. 
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Figure 5.5. Modal frequency response of the single layer 
configuration as computed in NASTRAN. 

[-: reference plate ,-: single layer] 
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The modal loss factors estimated using the modal strain energy 
method are compared to those measured experimentally for the single layer 
configuration in Figure 5.6. The estimated loss factors are greater than the 
experimentally determined loss factors throughout the spectrum of 
interest, and especially in the middle frequencies. 

The modal frequency response of the single layer configuration is 
compared to the experimentally measured frequency response in 
Figure 5.7. The comparison was accomplished by normalizing both the 
experimentally determined frequency response and the frequency response 
computed in NASTRAN. Both responses were normalized using a value of 
1.0 ^ n ^ C ' • The effects of the greater loss factors estimated by the modal 
strain energy method are obvious as the level of the predicted response is 
lower than the measured response. The shift in frequency between the two 
curves is due to the finite element model being inherently stiffer than the 
actual system. The correlation between the two curves is especially good 
below 250 Hz as this is where the differences between estimated and 
measured modal loss factors are the least. 
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Figure 5.6. Comparison of estimated and measured modal loss 
factors for the single layer configuration. 


« 

♦ a Q 


□ □ 


a • 


♦ 

□ 


a 15.6 °C 
♦ MOE Method 


QQ 


- 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 — 

0 200 400 600 800 

Frequency (Hz) 


5 9 





FREQUENCY RESPONSE (DB) 



Figure 5.7. Frequency response comparison for 
the single layer configuration. 
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C. DOUBLE DAMPING LAYER 


The modeling of the double constrained layer damping system was 
accomplished as shown in Figure 5.8. The double layer configuration 
consists of a base layer modeled with offset QUAD4 elements, two 
viscoelastic layers consisting of HEXA elements, and the top constraining 
layer modeled with offset QUAD4 elements. The middle constraining layer 
was modeled using three layers of HEXA elements in order to give this 
layer the stiffness necessary to act as a constraining layer. The model was 
meshed using 60 elements in each layer as shown in Figure 5.8. 

Using this model, the modal strain energy method was applied to 
determine approximate modal frequencies and loss factors for the double 
layer configuration. To determine the loss factors, normal mode 
extractions were performed using reference frequencies of 50, 200, 500, 800, 
and 1000 Hz. A composite set of modal loss factors for the double layer 
system was then compiled based on the estimated modal frequency’s 
relation to the reference frequency used to calculate modal strain energies. 
This composite set of modal loss factors is listed in Table 5.2, and is plotted 
in Figure 5.19 as a comparison to the experimentally determined modal 
loss factors for the double layer configuration. The estimated modal loss 
factors for the double layer show high damping, but they compare favorably 
with those measured experimentally. 









Modal Loss Factors 


TABLE 5.2. ESTIMATED MODAL LOSS FACTORS FOR THE DOUBLE 

LAYER CONFIGURATION. 
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Figure 5.9. Modal loss factors for the double layer configuration 
as determined from the modal strain energy method 
and determined experimentally. 
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The modal frequency response of the double layer configuration was 
computed in a manner similar to the single layer in that smoothed loss 
factor data was used in the MSC/NASTRAN damping table. Likewise, a 
unit excitation force was applied, and the first thirty modes were used in 
the modal summation. The results of the modal frequency response 
calculation are shown in Figure 5.10. The dashed line represents the 
undamped response, and the solid line represents the frequency response of 
the double layer configuration. 

The frequency response calculated using NASTRAN was compared 
to the experimentally determined frequency response of the double layer 
configuration as shown in Figure 5.11. Comparison of the two frequency 
response curves shows similarity in form, but a much lower response level 
for the numerically determined response. Once again this could be due to 
the higher damping predicted by the modal strain energy method and the 
inherently higher stiffness of the finite element model. 
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Figure 5.10. Modal frequency response for the double layer 
configuration as computed in NASTRAN. 

[-: reference plate ,-: double layer] 
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Figure 5.11. Comparison of the experimentally determined 
and numerically predicted frequency responses 
for the double layer configuration. 
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D. POCKET PLATE RESULTS 


The pocket plate configuration was modeled with the same offset 
plate elements and solid viscoelastic elements as the single layer 
configuration. However, the pocket plate required the modeling of the 
milled structure around the viscoelastic material and the welds between 
the cover plate and milled plate. A representation of the model is shown in 
Figure 5.12. The base structure, cover plate, and the structure immediately 
around the cover plate was modeled using offset plate elements. The 
viscoelastic material and the portion of the structure immediately adjacent 
to it were modeled using the solid HEXA elements. Since the viscoelastic 
material and cover plate are physically separated from the surrounding 
plate, except where the viscoelastic is adhered to the base structure, care 
was necessary in creating the finite element mesh. 

The model was created in PATRAN using PATRAN’s node editing 
and equivalencing capabilities [Ref. 18]. This allowed the generation of a 
finite element mesh with two nodes at the same geometric point in space. 
Using this node editing technique, a mesh was created which allowed the 
viscoelastic and cover plate to vibrate separately from the surrounding 
structure, yet at the same time, keep the number of elements and nodes in 
the model to a minimum. The welded points on the cover plate were also 
modeled using node editing techniques. At weld points the finite element 
node on the cover plate was equivalenced with its corresponding node on the 
base structure, resulting in a single node and a connection between an 
otherwise separate base structure and cover plate. At non-welded points 
on the cover plate there were two nodes at the same geometric point; one to 




represent the cover plate, and the other to represent the base structure. The 
model was meshed using a 5x11 mesh resulting in 40 elements in the cover 
plate, viscoelastic and base layer as shown in Figure 5.12. 

Using this modeling scheme, the modal strain energy method was 
employed to estimate the modal frequencies and loss factors. Normal mode 
extractions were made using viscoelastic material properties at 50, 200, 500, 
800, and 1000 Hz. Using these reference frequencies a composite set of 
modal loss factors was obtained. These loss factors are listed in Table 5.3 
and are plotted versus frequency in Figure 5.13. The estimated loss factors 
give good damping over the spectrum of interest with an average modal loss 
factor of 0.075. As with the previous cases, damping values used for the 
modal frequency response calculation came from a curve-fit to the set of 
composite modal loss factors. 

The modal frequency response of the pocket plate was computed 
using the first 30 modes and viscoelastic material properties at 200 Hz. The 
resulting estimated frequency response is shown in Figure 5.14. The 
response shows a definite frequency shift to the left along with good 
damping of the frequency response when compared to the undamped 
response. 
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Modal Loss Factor 


TABLE 5.3. ESTIMATED MODAL LOSS FACTORS FOR THE POCKET 


PLATE CONFIGURATION. 
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Figure 5.13. Estimated modal loss factors versus 
frequency for the pocket plate. 
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for the pocket plate configuration. 


(-: reference plate 


: pocket plate] 















The estimated modal loss factors and modal frequency response for 
the pocket plate were compared to those measured experimentally. The 
loss factor comparison is shown is Figure 5.15 and the frequency response 
comparison is shown in Figure 5.16. The estimated modal loss factors are 
higher than those measured experimentally, however, the frequency 
responses compare quite favorably with each other. The frequency 
response curve calculated through finite element analysis has a lower 
response level and a frequency shift to the right of the measured frequency 
response. This is expected due to the increase in damping predicted by the 
modal strain energy method and by the fact that the finite element model is 
inherently stiffer than the physical system. 
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Figure 5.15 Comparison of experimentally measured and numerically 
estimated modal loss factors for the pocket plate configuration. 
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Figure 5.16. Comparison of experimental and predicted frequency 
responses for the pocket plate configuration. 


73 







VL CONCLUSIONS 


Constrained viscoelastic layer damping is an extremely effective 
method for reducing broadband vibration. In each of the experimental 
cases the peak amplitudes of frequency response were reduced by 
approximately 25 decibels, a reduction of 18 times below the undamped 
reference plate. 

Of particular note are the performance of the pocket plate and 
floating element configurations. Although they are not ideal 
configurations in terms of “true” constrained viscoelastic layer damping, 
the damping levels achieved are quite satisfactory. In a final comparison of 
experimentally determined modal loss factors for the four treatments, the 
double layer configuration yields the greatest damping and the pocket plate 
the least. This comparison is shown in Figure 6.1. As was previously 
reported in Section IV, tl floating element configuration yielded an 
average increase of 25 percent over the modal loss factors of the pocket 
plate. It is also noted that the average modal loss factor for the pocket plate 
is approximately 50 percent of the average modal loss factor of the single 
layer configuration. Similarly, the average modal loss factor for the 
floating element configuration is approximately 50 percent that of the 
double layer treatment. 
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Figure 6.1. Experimental modal loss factor comparison 
for the four damping configurations. 


The modal strain energy method tends to overpredict the modal loss 
factors for these highly damped, thick plates. Although the modal strain 
energy method predicted modal loss factors greater than those measured, 
the same relative differences in modal loss factor between damping 
configurations are maintained as shown in Figure 6.2. The estimated 
modal loss factors for the pocket plate are approximately 60 percent less 
than those of the single layer configuration This indicates that there is a 
consistency between the damping values predicted by the modal strain 
energy method and those of the actual physical system. 
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Figure 6.2. Comparison of numerically estimated modal loss factors for 
the single layer, double layer, and pocket plate configurations. 
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There are several possible reasons for the differences between the 
experimentally determined modal loss factors and those estimated by the 
modal strain energy method. The first is that the material properties of 
ISD-112 as reported by the 3M Corporation on the reduced frequency 
nomogram may not be consistent with the material properties actually 
present in the material used. Since the numerical analysis was based on 
the reported material properties this is a possible source of uncertainty in 
the results. 

The second source of uncertainty is in the adhesion of the ISD-112 to 
the aluminum plates. Although the plates were clean when the ISD-112 
was applied, it was noted that the adhesive qualities of the ISD-112 were not 
uniform throughout the material. A lack of adhesion may cause a 
reduction in the damping capability of the system. 







Another source for the difference between the experimental results 
and numerical results lies within the finite element model. The number of 
elements used in the model greatly affects its “stiffness.” As the number of 
elements is increased the model should become less stiff and results are 
expected to approach those that are measured. Also, the type of element 
used to model the base layer and constraining layer may have an effect. In 
this research plate elements were used to model both the base and 
constraining layers. It is possible that one or more layers of solid elements 
may produce results that agree better with experimental results. 

The modal strain energy method is a useful tool for vibration 
damping design, however, the method appears to overestimate modal loss 
factors for thick, highly damped plates by as much as 50 percent. In each 
configuration that was modeled, the modal loss factors estimated below 
300 Hz were fairly close to the actual values, and above 300 Hz the difference 
between estimated as measured loss factors increased. 

One drawback to the modal strain energy method in design is the 
CPU time required for normal mode extraction and modal frequency 
response; especially in complex structures with a large number of 
elements. Although the modal strain energy method is good for finalizing 
a design, it should not be used for design optimization due to the amount of 
CPU time required for normal mode extraction and modal frequency 
response calculations. 
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VIL RECOMMENDATIONS 


The thick plate used in this research is a generic model of many 
physical systems that may see use in naval application. There are several 
areas which deserve more research and clarification, including the 
following: 

• Model the floating element configuration in finite elements and 
check the effectiveness of the modal strain energy method in 
predicting modal loss factors for this configuration. 

• Investigate the capability of the modal strain energy method to 
predict loss factors for thick, lightly damped plate configurations. It 
is possible that the modal strain energy method works better for 
lightly damped configurations. 

• Investigate methods for improving the adhesion of ISD-112 to the 
base structure and constraining layer. 

• Investigate the possible use of the floating element as an additional 
vibration absorber by changing its mass in order to damp a specific 
frequency. In this same idea, the use of tuned dampers along with 
the constrained layer treatments could be looked at. 
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APPENDIX A 


FORTRAN PROGRAM USED TO COMPUTE MODAL LOSS FACTORS 
FOR THE SINGLE DAMPING LAYER DESIGN 

This program used the Ross-Kerwin-Ungar equations of Section III 
to compute estimated modal frequencies and loss factors for the single layer 
configuration. Modal frequencies for an undamped plate are read from a 
data file and estimated modal frequencies and loss factors for various layer 
thicknesses are output to another file. Material properties of ISD-112 are 
computed using University of Dayton data and curve-fitting equations to the 
reduced frequency nomogram [Ref. 9,12]. The units used in this program 
are pounds, inches, and seconds. 
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PROGRAM CARTE! 

*xxxxxxx::xxxxxxxxxxxxyyxxxxxxyxxxxxxxxxxxxxyxx*xxxxxx*x*x*xxxxxxxxxxx 
THIS PROGRAM IT, TO CALCULAI E SYSTEM LOOP: FACTORS FOR VARIOUS 
T EHF ER A I URE S All 0 MOOES OF A COtlS I RAINED LAYER VISCOELASTIC 
DAMPING SYSTEM USING THE ROSS-KERl11II-UIIGAR EQUATIONS. SYSTEM LOSS 
FACIORS ARE COMPUIEO TOR A IEMP E RA TORE RANGE OF 50-100 DEGREES 
FARENNLIT. BASE TLAT E THICKNESSES VARY FROM 0.575 TO 0.75 INCHES 
III 0.125 1NCREMEII! S . FOR EACH BASE PLATE THICKNESS, THE VISCOFLASTIC 
THICKNESS IS VARIED TROM 0.015 TO 0.060 INCHES 111 0.015 INCH 
INCREMENT S. FOR EACH VISCOELASTIC THICKNESS THE CONSTRAINING 
LAYER THICKNESS IS VARIED FROM 0.0625 TO 0.25 INCHES IN 
0 0625 INCH INCREMENTS. 

THIS PROGRAM AFLIES TO A FREE-FREE-FREE-TREE PLATE. 

THE VISCOELASTIC MATERIAL IS 5M ISD-112. 

THE FULL 01II III! COEFFICIENTS ARC DEFINED: 

El - YOUNG'S MODULUS OF BASE PLATE (PSI) 

E2 = YOUNG'S MODULUS OF VISCOELASTIC LAYER (PSI) 

E3 = YOUNG'S MODULUS OF CONSTRAINING LAYER (PSD 
02 - SHEAR MODULUS OF VISCOELASTIC MATERIAL (PSI) 

IIUI = POISSON'S RATIO OF BASE PLATE 
IILi2 - POISSON'S RATIO OF VISCOELASTIC MATERIAL 
111 - THICKNESS OF BASE PLATE (IN) 

IIP = I HI CKHESS OF VISCOELASTIC LAYER (111) 

II 5 = THICKNESS OF CONSTRAINING LAYER (IN) 

HIOT = TOTAL PLATE THICKNESS (IN) 

R11U1 = DENSITY OF BASE ELATE (LBF-SECXX2/INxxQ) 

RU02 = DENSITY OF VISC0ELAS1JC MATERIAL 
RIIOS = DENSITY OF CONSTRAINING LAYER 

I 0 , T ROM, ITROM . II, ML , ET FROL , SL , SH , TROI. , 8 C ARE COEFFICIENTS 
FOR THE REDUCED FREQUENCY NOMOGRAM EQUATIONS 
FP - MODAL FREQUENCY OF THE UNDAMrED BASE PLATE (HZ) 

TCP = MODAL FREQUENCY OF THE COMPOSITE PLATE (HZ) 

EI A 2 = LOSS FACTOR OF VISCOELASTIC MATERIAL 
EI AS = SYSTEM LOSS FACTOR 
KQR - HAVE NUMBER OF PLATE 
GC = GRAVITATIONAL CONSTANT (IN/SECxxP) 

I = TEMPERATURE OT VISCOELASTIC MATERIAL (DEG F) 

TIP = MODAL FREQUENCY OF UNDAMPED PLATE (RAD/SEC) 

TP - MODAL FREQUENCY OF UNDAMPED PLATE. (HZ) 

HIT = MODAL FREQUENCY OF DAMPED PIAIE (RAD/SEC) 

TCP e riODAL TREQUENCY OF DAMPED PLATE (HZ) 

THE UNITS USED IN THIS PROGRAM ARE LB, INCH, SEC, DEGREES F 

UIIDAfllTD MODAL FREQUENCIES ARE FROM A FINITE ELEMENT ANALYSIS 
OT THE I REE-FP.EE-FREE-FREE TLATE. 

RESULTS APE OUTPUT TO DATA FILE "LOSFCTR DATA" 

PRIOR TO RUNNING THE PROGRAM TYPE THE COMMAND 
"IILEDEE LOSrCTR DISK LOSTCTR DATA” 

FINITE EIEMCNT MODAL PREQUENCIES ARE INPUT FROM TILE: 

"FL TrRQ DATA” 

XXXVXX*X*XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXyXXXXXXXXXXXXXXXXXX 

real El . E2, E3, NUI , NU2, HI , H2, H3, RIIOl . RH02 , RH03, NT 0 T 
REAL ID, [ ROM,MROM,N,ML,EIrROL,SI,SH,FROL,C 
REA I TP,FCP,ETA2.ETAS,KQR,PI,GC,T 
REAI FR10,ETA21D,SUP1,SUB2,SUB3,M10 
REAL HP 1 , H 7 ,l , G. Cl ,ALPHRE. ALrillM, BRE, DIM, DEL RE 
REAL DEI 111, EHCUBE, 5UB6, SUB5, HP,HCP, DENS, SUB6 , SUB7 , SUBS 
IIII EGER V 
DIMENSION TIP(S) 

PI = •'*. XA1 AIK 1 . ) 

OrCN(IINII-10,rilE-'LOSFClR',5TATUS-'OID') 

OPETK UNIT =11,n L L='PLITRQ',SIA1US^'OLD') 

C ASSIGN MATERIAL CONSTANTS FOR BASE TLATF,CONSTRAINING LAYER 
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C Ann VISCOELASTIC 

c 

F 1 = 1 . 0 E 7 
HI=0.750 
RHO 1 = . 0168 
lllll = . 33 
RHOc-.035 
E 3 =1 .0E7 
RHO 3 -.0968 
OC=386. 

IHJ 2= . 5 
C 

C DCriNE COMSTAIITS FOR NOMOGRAM EQUATIONS 

C 

T 0 = 1 0 6 . 0 
FROM-2.OF6 
HROM-688.96 
II-.275 
ML =3.7 
El FROL =1 . 08 
51=0.65 
5H = -.55 
FROL = 5000 . 0 
C = 2.5 
C 

C BASE FLATE THICKNESS LOOP 
C 

111 = 0.375 
DO 100 J = 1,6 

HRITE(10,703) 'BASE PLATE THICKNESS, H1=',H1 
703 F ORMAT(A27,F6.3) 

MR 1rc<10,x) 

C 

C VISCOELASTIC THICKNESS LOOP 
C 

(12 = 0.015 
DO 150 V = 1,6 
C 

C f 0113 1 RAI I/I f/G LAYER THICKNESS LOOP 
C 

H3 = 0.0625 
DO 20 0 L = 1 , 6 

IIRITE( 1 0,7 02 ) 'VISCOELASTIC LAYER THICKNESS, H2=’,1I2 
UR I I EC 10,XT 

HR I TEC 1 0,7 02 ) 'CONSTRAINING LAYER THICKNESS, H3=’,H3 
7 02 FORMA 1(A35,F6.6) 

HR I 1C(10,X) 

HR 11£1 10,700) 'TEMP','MODE','FCP','ET A2','ETAS', 'G2' 

700 FURMAT(A3,3X,A5,6A15) 

C 

C PLATE MODE LOOP, READ MODAL FREQUENCY AND COMPUTE HAVE NUMBER 

C 

HT 0T = H! HI2HI3 
DO 300 1=1,8 
RE ADC 1 1 , x ) lirci) 

C 

C CALCULATE MODAL FREQUENCY AND HAVE NUMBER OF BASE PLATE 
C 

SUB 7 = SQRT((El*(HTOTxx3)xGC>/<12.*C1 .-IIU1 XX2 ) XRH01*H1)) 
KQR-HPtIJ/SUB7 

r p=hp( i)/( 2 .xrn 
c 

C TEMPERATURE LOOP 
C 

T = 30.0 

DO 600 K=1,15 
C 

C CALCULATE PROPERTIES OT VISCOELASTIC FOR GIVEN TEMPERATURE AND 
C MODE 
C 

501 FR10 = l 0010(PP)-U2.MT-T0))/( 525.41-10) 

FR = 10. «*(FR 10) 
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A-oRio-iooiocrRom/c 

SUB1-C*((SL* SH>xA K SL-SH)*(1.-5QRH1.+AXX2.)))/2. 
El A21U=LOG10( El FROL ) fSUBl 


VISCOELASTIC LOSS FACTOR 

EIA2=10.xx(ETA210) 
SUB2=2.XLOG10(MROM/HL) 

5 U B 3 = 1 . » ( FROM/ FR ) x XN 
Ml 0 = L 0G10(ML1+3UB2/SUB3 

VISCOELASTIC SHEAR MODULUS 


G2=10.*x(M10) 


VISCOELASTIC YOUIIG'S MODULUS 


E2-G2*2.x(l .h MU2) 

CALCUI AT IOIIS FOR DAMPED PLATE AT TEMPERATURE T USING ROS 
KERCH M-UNGAR EQUATIONS 


0-G 2/(E5XH3XH2XKQR) 

U21 = (HHH2)/2. 

HS1=H2* (111 +M3 )/2 . 

Cl - Cl*lll*( 1. 4G) »GXE3XH3 
D-G>ETA2X(E1XH1+E3XH3) 

AI PURE = G*E1 xill XE3XH3X( 1131 xx2 . )x( Cl +D*ETA2) 

ALrill 11-GxE1XH1XE3XH3X11131XX2. )*C1*ETA2 
PRE = E1XII1XE2*H2*H31X(C1 + DXETA2) 
B1I1=E1*II1XE2XH2*H31X(C1 XETA2-D) 

SUB8 = 2.XG*E2*H2*E3XH3XU21 *H31 XC 
DELRE=SIJB8x((l.-(ErA2xx2.))+(DX(2.xfTA2))) 
DFLIM-SUBSx(Clx(2.XEIA2)-(Dx(l.-(tTA2X)f2.)))) 
:>UBA = ( 12 ./(Cl**2 . i Dx»2. 1 1 *(AlPHRE-BRE-DEERE) 
EHCUBE r (El*(Hl*x3.))*E3*(H3**3 . MSUBQ 

MODAL FREQUEIIC f OF DAMPED PLATE 

DENS = RH01 x|ll-*RH02»H2i HJ*RM03 
SUB5=(EHCUBEXGC)/(12.xn.-HU1XX2.1XDEMS) 
ITCP-F.QRxSQR i ( SUB5 ) 
rcP = HCIV(2.>P1 ) 

COMPARISON OF FT AND TCP 

ir (abs( i .-rr/rcp) .le. o.n then 

GO I 0 500 
ELSE 

FP-rcp 

GOTO 501 
END1F 

COMPUIE SYSTEM LOSS FACTOR 

00 5UB6 = (12./(Cl**2.-tn* v 2. )) x ( AL PMIM-BIM-DELIM) 

E I AS - ( 1 . / EHCLIBE ) XSUB6 

PRINT RESULTS 

HR IT E(10,7 01 ) T , I ,FCT,ETA2,ETAS,G2 
01 FORMAT(5X,F7.3,2X,I 2,3X,0E15.A ) 

NEXT TEMPERATURE 

T - T ♦ 5.0 
00 CONTINUE 

NEXT MODE 


HR I T E(10,* 1 
300 CONTINUE 







o n >— nnnro o o n 


RCl-IINDC UNI T = 11 > 


NEXT CONSTRAINING LAYER THICKNESS 

H 3 = H 3'0.0625 
00 CONTINUE 

NEXT VISCOELASTIC LAYER THICKNESS 

N2 r N2'0.015 
50 CONI INUE 

NEXT BASE rLAT E THICKNESS 


Hl^Hl+n.125 
100 CONTINUE 

ClOSE(NNIT'IO) 
ClOSEC UNIT = 11) 
END 
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APPENDIX B 


FORTRAN PROGRAM USED TO COMPUTE MODAL LOSS FACTORS 
FOR THE DOUBLE DAMPING LAYER CONFIGURATION DESIGN 

This program uses the Ross-Kerwin-Ungar equations of Chapter 3 
to compute estimated modal frequencies and loss factors for the double 
layer configuration. Modal frequencies for an undamped plate are read 
from a data file and estimated modal frequencies and loss factors for 
various layer thickness combinations in the double layer configuration are 
output to another data file. Material properties of ISD-112 are computed 
using University of Dayton data and curve-fitting equations to the reduced 
frequency nomogram [Ref. 9,12]. The units used in this program are 
pounds, inches, and seconds. 
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no o o oooooonooooooon ooooooooooooonooonooooooooonnoooooonoooononn 


PROGRAM 1 MOLYR 


XXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

THIS PROGRAM COMPUTE': THE SYSTEM I OSS FACTOR AND MODAL FREQUENCY OF 
A DOODLE CONS I RAJ NED LAYER V ISCOELASTICALLY DAMPED PLATE. THE LOSS 
FACTORS ARE COMPUTED FOR A SPECIFIC PL A IE/DAMP ING LAYER 
CONFIGURATION AND OVER A TEMPERATURE RANGE OT 30-100 DEGREES 
FARENIIE1 T . 

THIS PROGRAM APPLIES TO A FRFF-FREE-FREE-FREE PLATE AND THE UNDAMPED 
MODAL FREQUENCIES ARE DETERMINED FROM FINITE ELEMENT ANALYSIS AND 
ARE READ INTO THIS PROGRAM FROM FILE 'PLIFRQ*. 

THE VISCOELASTIC MATERIAL IS 3M ISD-112. VISCOELASTIC MATERIAL 
DATA IS FROM UNIVERSITY OF DAYTON RESEARCH INSTITUTE 

THE UNITS USED III THIS PROGRAM ARE ID, INCH, SEC, AND DEGREES F. 

LOSS FACTORS, DAMPED FLATE MODAL FREQUENCIES, AND ISD-112 FROTERTIES 
ARE COMPUTED III SUBROUTINE 1 R K U' FOR EACH ll-TH CONSTRAINED LAYER 
SYSTEM. UPON COMPLETION OF EACH TEMPERATURE COMPUTATION THE SYSTEM 
LOSS FACTOR AND CORRESPONDING DAMPED TLATE FREQUENCY ARE MR IT T Ell TO 
FILE 'T MOLYR DATA' . 

THE F 01. L ON IIIG COEFFICIENTS ARE DEFINED: 

EL = YOUNG'S MODULUS OF HASP PL A1E (PSI) 

E3 = EQUIVAIENT YOUNG'S I10DUIUS OF N-TH CONSTRAINED LAYER SYSTEM 
E3PP.M = YOUNG'S MODULUS OF IHF CONSTRAINING LAYER IN N-TH LAYER 
EIICURE = EQUIVALENT STirrilESS OF CONS I RAIHED LAYER SYSTEM 
ETAS - SYSTEM LOSS EACIOR 

E T A j - LOSS FACTOR OT NTH CONSTRAINED LAYER SYSTFM 
EIA 3 FIT -- LOSS FACTOR OF CONSTRAINING LAYER IN THE N-TH LAYER ( = 0 
I CP - FREQUENCY OF IIIE DAMPED PLATE (HZ) 

FP = FREQUENCY Of UNDAMPED PLAIE (HZ) 

GC = GRAVITATIONAL CONSTANT 

111 - THICKNESS OF BASE FLAIE 

112 = THICKNESS OF 1ST VISCOELASTIC I AYER 

Hi = 1HICKIIESS OF ll-TH COIIS I R A1 HE D LAYER SYSTEM 

HA = T 0 T Al DAMPED PLATE THICKNESS 

HITRIT = THICKNESS OF BASE LAYER III N-TH LAYER 

H2PRM - THICKNESS OF VEIT IN N-TH LAYER 

H3PRM - THICKNESS OF CONS 1 RAI III IIG LAYER IN ll-TH LAYER 

KQR = HAVE tIUMDCR OF UNDAMPED r I ATE 

H "1 = roiSSOII'S RATIO OF BASE FIATF AND CONSTRAINING LAYERS 
1102 - roiSSOII'S RAIIO or VISCOELASI IC MATERIAL 
HU2PRM = POISSON'S RATIO OF VEIT IN ll-TH LAYER 
R1101 = DENSITY OT BASE Ft ATE AND CONSTRAINING LAYERS 
RIIU2 - DENSITY OT VISCOELASI IC MAIERIAL 

Riiuzrn = ocNsiir or vem in n-ih layfr 

RII03 = DENSITY OF N-TH LAYER ( = RH01) 

RII03FM = DENSITY OF ll-TH LAYER COIIS IRAI HI NG LAYER 
T - TEMPERAIURF VARIABIE 

MP = FREQUENCY OF UNDAMPED TLATE (RAD/SEC) 

C.F.IFROI . r ROL , F ROM . ML , MROM, N , SH, SI , T 0 = COEFFICIENTS FOR THE 
REDUCED FREQUENCY NOMOGRAM EQUATIONS 

XXX*X**VXX**XXXXXXX*X*XXX¥XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

REAI C, FI . E 3 , E 3 T R. M, EHCUBE, El AS , E T A 3, V IA3PM, ET T ROL , TCP, FP, FROL 
REA I TROM.GC,HI,H2,M3.HA,HirRM, II2PRM,H3PRM.KQR,ML,MROM,N,NU1 
REA1 NN2 , IIU 2 PM, r I , RHUi , RII02,RH02rM,SH,3L , SUB1 , T , T 0 , IIP 
REAL P.H03, P.H03PM 

DIMENSION Hr(17) 
ri-A.OxAIAN(l.O) 

OPLII(HIMT-]0,ril E-'PI irRQ' ,STATUS-'01 D' ) 

OPEIK UNII-ll,riLE-'IHULYR',STATUS-'OLD' ) 

ASSIGN LAYER THICKNESSES 
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n n n non non n r> n -j o o o n n n o non 


111 -0 . 5 

112 = 0 . 065 
HH'RH=0.09375 
M2f'RH = 0 .015 
H3PRM=0 . 09575 

115 = 111 I'RM 1112PKM * H5PRM 
Mm = It 1 i H2 < 113 

A55IGII MATERIAL C OH 3 I All 7 S 

E 1 = 1 . 0 E 7 
RIIU1 =0.0968 
mil =0.33 
RH02- 0 . 035 
11112 = 0 . 59 
RII02PM = 0.035 
m)2r’M=0 . 5 
RH03-0.0968 
E 5 PR M=1.0E7 
RII0;PM = 0 . 0968 
ETA 3PM 1 0.0 
El A3 = 0.0 

GC- 386 . 0 

DEEIHE VISCOELASTIC COMSTAMTS TOR NOMOGRAM EQUATIONS 

10= 106.0 
T ROM-2.0E6 
i:R(.'M = 638.96 
11 = 0.27 5 
ML =8 . 7 
E I f ROL-1.08 
SI -0.65 
311-0. 55 
F RUl =5000.0 
C = 2.5 


WRITE TLAIE CHARACTERISTICS 


00 
0 1 


IIRI im 1 ,700 ) 

hrii n11,700) 

IIRI If ( 1 1 . 7 00 ) 
IIRI I r.U I . 7 00 ) 
I IK I I r C J 1 ,700) 

itri i e< 11 , n 

muim i, x) 


111 r ’ , H1 
(12 = • .112 
II1 PRI ITE = 
II2TRI ITE = 
IISTR1 ME = 


II1PRM 
II2IRM 
II SERfl 


I ORMAI(A 12,E6.5 ) 

IIRI 11 C 1 I , 701 ) ' ( EMP ' . 'MODE' , TCP 
FORMA I(AS,3X,A5.2A15) 

hr ifEC1l,* ) 


'ETAS' 


PLATE MODE LOOT 


DO 100 1-1,17 
READ! 10, X) IIP(1 ) 


CALCU1AIE MODAL FREQUENCY AND WAVE NUMBER OF UNDAMPED PLATE 

Sum - SORT ( ( El*(H6x>*3)*GC)/U2 .*( 1 . - MU 1 x *2 )XRHO1 *116 ) ) 

K Q R -11P ( I ) / S U B1 
r P = lll ( I )/( 2 TPI ) 


T EIIPERAI URE LOOP 
I - 50 . 0 

DO 20 0 KT . 15 

COMPUIE SYS I EM LOSS TACTOR AND STirFNESS FOR N-TH LAYER 

CAL L RKUC 111 PRM, II2PRM, H3PRM, El , RHO1 , HU 1 , NU2PM, RHOZHT, E3TRM, RH03PM, 
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cri a sen, i, rr, kqr, etas, eiicube, rcr, to,proh.mrom,ii,ml, ETFROL,3L,SH, 

C P R 01 , C ) 

c 

C CONYERI RESUL T 3 PROM II - TII LAYER CALCULAriOIl 10 TOTAL PLATE 

c 

E T A 3 - 0.0 

FP = UPC I)/(2.*PI ) 
c FP-rr.p 

C 3 = ENCODE/ ( H3*x3 ) 

C 

C COMPUTE SYSTEM LOSS PACT OR AMD FREQUENCY FOR TOTAL PLATE 
C 

CALI. RKtK 111 ,112, U3, El . RIIOl , Mill , IIU2, RH02, E3, RH03, ET A3, T, FP, KQR, ET AS, 

ccncuBE,r cp,t o,rRon, drum,n,ml,eifrol,sl,sii,prol,C) 

c 

C rRIMT RESULTS 

URITEL11,702) T,l,PCr,ETAS 
7 02 FORMAT(5X,F7.3,2X,I 2,3X,2E15 . A ) 

C 

C NEXT TEMPERATURE 
C 

T - T i 5 . 0 

roo con ri hue 

c 

C NEXT MOPE 
C 

HR IT E(11,X ) 

100 CONTINUE 
C 

CLOSECIINI 1=11) 

Cl USE(UNI 1=10) 

END 

C 

o «**'»*xxx*vxxxxx*x»xxy»vxxxxxxxx*xx'«'*xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C XXXXXXXXoXXyXXKXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

SUBROUTINE RKU( HI , H? , M3, E1, RHOl • Mill, IIU2, RII02, E3 , RMO 3, E T A 3, T , FP , 

C KQR.ETAiS,El (CUBE,TCP,TO,PROM, MROM,II, ML,ETFROL,SL,SH,FROL,C) 

C 

C XX*VXXXXXXX*XXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

C THIS SUBROUTINE CALCULATES VISCOPIASTIC PROPERTIES BASED ON THE 
t: UN I VERS IT ( OP DAYTON REDUCED FREQUENCY NOMOGRAM EQUATIONS, AND THEN 

C CALCULATES PLATE STIEPNESSES AND LOSS FACTORS BASED ON THE ROSS- 

C KERIT1 N-UNOAR EQUATIONS. 

C 

C THE rnunilHIG ADDITIONAL variables are DEFINED TOR USE in this 

C SUBROUTINE: 
c 

c a = cnrrnciENT tor nomogram equations 

C ALPIIIM,ALPHRE = IMAGINARY AND REAL COMPONENTS OF COEFFICIENT 

C ALPHA IN THE KKU EQUA UNITS 

C Dill, BP,E = IMAGINARY AND REAL COMPONENTS OF COEFFICIENT 'B' IN 

C THE RKU EQUATIONS 

C Cl.D - COET riCIEIUS FOP. RKU EQUATIONS 

C DELIM,DELRE = IMAGINARY AND REAL COMPONENTS OF COEFFICIENT 

C DELTA IN THE RKU EQUAIIOTIS 

C DENS = cnnniNAlION or MATERIA! DENSITIES U5ED TO COMPUTE THE 

C FREQUENCY OE THE DAMPED PLATE AND DAMPING LAYERS 

C EIICUBE = EQUIVALENT STIEENESS OF DAMPED PLATE AS COMPUTED 

C USING RKU EQUATIONS 

C FTA2 = LOSS FACTOR OE V EM COMPUTED IN NOMOGRAM EQUATIONS 

C E I A210 - LOGIOTFIA?) 

C ETAS - SYSTEM LOSS FACTOR COMPUTED BY RKU EQUATIONS 

C FCT - MODAL PRCOUFNCY OP DAMP C D PLATE ( 117 ) 

C FP = MODAL FREQUENCY OP UNDAMPED PL A IE (HZ) 

C I R = REDUCED PREQUEIICY OP VEM 

C PR10 - L OG1 ()( PR ' 

C G2 - SHEAR MODULUS OP VEM 

C HZ 1 . N '> 1 = RKU EQUATION COErElCIFNTS 

C nin = L OGl 0 (G2) AS COMPUTED BY NOMOGRAM EQUAIIOTIS 


87 






non non n non non onnoon mono o noon 


lice - FREQUENCY or DAMPED PLATE (RAD/SEC) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

REAL A.AI PHIM,ALPURE,DIM,DRC.C.Cl,D,DELHI,DEI RE,DENS,El,E2,E3 
RLAl EHCUBE, FT A?, El A3, ET L21 0, E 1 AS, E J FROL , TCP, FP , TP., FRIO, FROL , FROM 
REAL G.G2,PC,H1 ,112,113,1121 , H3 1 , KQR, II1 0 , ML , HKOM, II, IIU1 , IIU2, PI ,RH0] 
RLAl RII02, RIIO 3, 311, SI , SUB 1 , SUD2, SUB 3, SUB A, SUB5,5UB6,SUB 7 , SUD8 
RLAL SUD9.SUU10,3UB11,T,T 0 ,IICP 

PI-A.OXATAIK 1.0) 

GC= 386.0 

CALCULATE PROPERTIES OF VEM TOR GIVEN TEMPERATURE AMD MODE 

1 FR10=LOG10(FP)-(12.x(T-T0))/(525.-H-T0) 

FR"10.xx(tR10) 

A = ( | R1 0 — u(JOlot FROL ) )/C 

SUDl=Cx((SL+SH)xA+(SL-SH)x(l.-SORT(1.+A*x2)))/2. 

ETA210-LOG10(ETFROL)+SUDl 

VISCOELASTIC LOSS FACTOR 

E T A 2 = 1 0 .xx(ETA210) 

VISCOELASTIC SHEAR MODULUS 

SUB2 = 2 . xL0G1 0(MROIVML ) 

SUB3 = 1 . + C FROM/FR)xXU 
M10 r LOG10(ML)+SUD2/SUB3 
G2 = 10.*x(Ml 0) 

VISCOELASTIC YOUNG'S MODULUS 
E2 r G2*2.*(l .-MIU2) 

CALCULATIONS FOR DAMPED PLATE USING RKU EQUATIONS 

G-G2/( E3XH3XH2XILQR) 

HZl=CHl+H2)/2. 

M31 = H2+(HI + H3)/2. 

Cl=El*Mlx(l.+G)+G*E3xH3x(1.-ETA2XETA3) 

D = GxElxt|] XETA2<G*E3*H3x( ETA2+ETA3) 

SUDA = G*EUHlxE3xH3x(H31*x2) 

ALPHRE = SUBAX(C1X(1.-ETA2xETA 3) + Dx(ETA2+ETA3) ) 

ALPHHHS'TBAxcci x( ETA2TETA3)) 

SUB5-El*HlxE2*H2XH31 

BRE-SUB5X(C1+D*ETA2) 

BIM=SUB5X(C1XETA2-D) 

SUD6=2.XGXE2XH2XE3*H3XH21XH31 
SUB7=1.-2.XETA2*ETA3-(ETA2XX2) 

SUB8=2.XETA2+ETA3-C ETA2xx2)xETA 3 
DELRE-SUB6x(ClxSUB7 + DxSUB8 ) 

DELIM-SUB6XCC1XSUB8-DXSUB7) 

SUBR=(12./(C1«x24Dxx 2))x(ALPHRE-BRE-DELRE) 
EHCUBE=Elx(Hlxx3)tE3x(H3xx3)+SUB9 

MODAL FREQUENCY OF DAMPED TLATE 

DENS-RM01XM14RH02XH24RH03XH3 
SUBI0=(EHCUBExGC)'(12.x(NU1XX2)xDENS) 

HCP r KQRxSQRT(SUOIO) 

F CP r IICP/ (2 . xri ) 

COMPARISON OF FP AND TCT 

Irc absc i.-FP/rcr) .le. o.io) then 

GOTO 500 

El SE 

Fp=rcp 

GOTO 501 
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nonn 


lifIDI F 


conruiE SYSTEM LOSS fACTOR 

SUBll-(12./(Clx*2*D*x2))*(AtPHIM“BIM-DElIM) 
ETAS=( 1 ./EHCUBE)*( (E3*ll3xETA3) + SUBll) 

C 

END 
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APPENDIX C 


DESIGN DRAWINGS FOR THE MACHINING 
OF THE FLOATING ELEMENT AND 
POCKET PLATE CONFIGURATIONS 


The drawings shown in Figures C.l and C.2 were used to machine 
the pocket plate and floating element plate used in the experiments. These 
drawings are included to show the relation of the pocket for the ISD—112 to 
the cover plate and how the viscoelastic was protected from the heat of 
welding. 
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Figure C.2. Design drawing of the floating element configuration. 
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APPENDIX D 


REPRESENTATIVE MSC/NASTRAN DATA DECK FOR 
THE DAMPING CONFIGURATIONS 

This data deck was used to compute the modal frequency response of 
the single layer damping configuration and is a representative sample of 
the NASTRAN decks used for the other finite element models. The values 
in the damping table are from a curve fit to the modal loss factors estimated 
from the modal strain energy method. Since the data deck for the normal 
mode and modal strain energy extraction is virtually identical to this deck, 
the Case Control deck commands for the normal mode analysis are 
included, but are commented out. 

The OUTPUT request provides data for an x-y plot of the modal 
frequency response. 

The units used in this deck are pounds, inches, and seconds. 
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Id single,■fr 
sol 30 
$ 

$ THIS DECK IS TO COMPUTE THE HODAI, FREQUENCY RESPONSE OF A 15i45 

$ INCH ALUMINUM PLATE WITH A CONSTRAINED VISCOELASTIC DAHP1NG LAYER. 

$ THE MODEL HAS 252 ELEMENTS WITH 84 ELEMENTS IN EACH LAYER. THE QUAD 
$ ELEMENTS ARE OFFSET FROM THE HEX ELEMENTS AS SUGGESTED BY THE 
$ LITERATURE. 

$ 

$ THE MATERIAL PROPERTIES OF ISD-112 ARE FROH THE 3H CORPORATION. 

$ 

$ THE MODAL LOSS FACTORS IN THE DAMPING TABLO ARE FROM A CURVE-FIT 

$ TO THE SET OF MODAL LOSS FACTORS COMPUTED FROH THE MODAL STRAIN ENERGY 
$ METHOD. 

$ 

$ THE UNITS USED IN THIS DECK ARE POUNDS,INCHES, AND SECONDS 

J I tiltllltlltltltlltlittllttltlltttltttltttlilttllllltlittllttttltltlilt 

tl«e 3000 
cend 

title = MODAL FREQUENCY RESPONSE / 200 HZ / 3M 

■ethod = 1 

spc = 1 

dload - 10 

frequency = 10 

sdasplng = 101 

set 111 - 95 , 

svector = all 

acceleratlon(plot,phase) = 111 
output{xyplot) 
xyprlnt acce / 95(t3) 

$ 

J litltttlttttltt 

$ THE FOLLOWING LINES ARE THE CASE CONTROL DECK CARDS FOR THE 
$ NORMAL MODE EXTRACTION AND STRAIN ENERGY REQUEST 
$ iethod = 1 
$ spc - 1 
$ set 10 = all 
$ set 11 = 169,thru,252 
$ ese = 11 
$ 

$ iiiiiiiniiiiti 

BEGIN BULK 


$ TITLE = SINGLE LAYER 

WITH QUAD OFFSET 

$ DATA DECK 

PRODUCED BY 

PATNAS VERSION 2.0: 24-NOV 

GRID 

1 

45.0000 15.0000 0.51500 

GRID 

2 

45.0000 15.0000 0.50000 

GRID 

3 

45.0000 12.5000 0.51500 

GRID 

4 

45.0000 12.5000 0.50000 

GRID 

5 

45.0000 10.0000 0.51500 

GRID 

6 

45.0000 10.0000 0.50000 


94 





GRID 

7 

45.0000 

7.50000 

0.51500 

GRID 

8 

45.0000 

7.50000 

0.50000 

GRID 

9 

45.0000 

5.00000 

0.51500 

GRID 

10 

45.0000 

5.00000 

0.50000 

GRID 

11 

45.0000 

2.50000 

0.51500 

GRID 

12 

45.0000 

2.50000 

0.50000 

GRID 

13 

45.0000 

0 . 

0.51500 

GRID 

14 

45.0000 

0 . 

0.50000 

GRID 

IS 

41.7857 

15.0000 

0.51500 

GRID 

16 

41.7857 

15.0000 

0.50000 

GRID 

17 

41.7857 

12.5000 

0.51500 

GRID 

18 

41.7857 

12.5000 

0.59000 

GRID 

19 

41.7857 

10.0000 

0.51500 

GRID 

20 

41.7857 

10.0000 

0.50000 

GRID 

21 

41.7857 

7.50000 

0.51500 

GRID 

22 

41.7857 

7.50000 

0.50000 

GRID 

23 

41.7857 

5.00000 

0.51500 

GRID 

24 

41.7B57 

5.00000 

0.50000 

GRID 

25 

41.7857 

2.50000 

0.51500 

GRID 

26 

41.7857 

2.50000 

0.50000 

GRIP 

27 

41.7857 

0 . 

0.51500 

GRID 

28 

41.7857 

0 . 

0.50000 

GRID 

29 

38.5714 

15.0000 

0.51500 

GRID 

30 

38.5714 

15.0000 

0.50000 

GRID 

31 

38.5714 

12.5000 

0.51500 

GRID 

32 

38.5714 

12.5000 

0.50000 

GRID 

33 . 

38.5714 

10.0000 

0.51500 

GRID 

34 

38.5714 

10.0000 

0.50000 

GRID 

35 

38.5714 

7.50000 

0.51500 

GRID 

36 

36.5714 

7.50000 

0.50000 

GRID 

37 

38.5714 

5.00000 

0.51500 

GRID 

38 

38.5714 

5.00000 

0.50000 

GRID 

39 

38.5714 

2.50000 

0.51500 

GRID 

40 

38.5714 

2.50000 

0.50000 

GRID 

41 

36.5714 

0 . 

0.51500 

GRID 

42 

38.5714 

0 . 

0.50000 

GRID 

43 

35.3571 

15.0000 

0.51500 

GRID 

44 

35.3571 

15.0000 

0.50000 

GRID 

45 

35.3571 

12.5000 

0.51500 

GRID 

46 

35.3571 

12.5000 

0.50000 

GRID 

47 

35.3571 

10.0000 

0.51500 

GRID 

48 

35.3571 

10.0000 

0.50000 

GRID 

49 

35.3571 

7.50000 

0.51500 

GRID 

50 

35.3571 

7.50000 

0.50000 

GRID 

51 

35.3571 

5.00000 

0.51500 

GRID 

52 

35.3571 

5.00000 

0.50000 

GRID 

53 

35.3571 

2.50000 

0.51500 

GRID 

54 

35.3571 

2.50000 

0.50000 

GRID 

55 

35.3571 

0 . 

0.51500 

GRID 

56 

35.3571 

0 . 

C.5COOO 

GRID 

57 

32.i429 

15.0000 

0.51500 
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GRID 

be 

32.1429 

15.0000 

0.50000 

GRID 

59 

32.1429 

12.5000 

0.51500 

GRID 

60 

32.1429 

12.5000 

0.50000 

GRID 

61 

32.1429 

10.0000 

0.5)500 

GRID 

62 

32.1429 

10.0000 

0.50000 

GRID 

63 

32.1429 

7.50000 

0.51500 

GRID 

61 

32.1429 

7.50000 

0.50000 

GRID 

65 

32.1429 

5.00000 

0.51500 

GRID 

66 

32.1429 

5.00000 

0.50000 

GRID 

67 

32.1429 

2.50000 

0.51500 

GRID 

68 

32.1429 

2.50000 

0.50000 

GRID 

69 

32.1429 

0 . 

0.51500 

GRID 

70 

32.1429 

0 . 

0.50000 

GRID 

71 

28.9286 

15.0000 

0.51500 

GRID 

72 

28.9286 

15.0000 

0.50000 

GRID 

73 

28.9286 

12.5000 

0.51500 

GRID 

74 

28.9286 

12.5000 

0.50000 

GRID 

75 

28.9286 

10.0000 

0.5i500 

GRID 

76 

28.9286 

10.0000 

0.50000 

GRID 

77 

28.9286 

7.50000 

0.51500 

GRID 

78 

28.9286 

7.50000 

0.50000 

GRID 

79 

28.9286 

5.00000 

0.51500 

GRID 

80 

28.9286 

5.00000 

0.50000 

GRID 

01 

28.9286 

2.50000 

0.51500 

GRID 

82 

28.9286 

2.50000 

0.50000 

GRID 

03 

28.9286 

0 . 

0.51500 

GRID 

84 

28.9286 

0 . 

0.50000 

GRID 

85 

25.7143 

15.0000 

0.51500 

GRID 

86 

25.7143 

15.0000 

0.50000 

GRID 

87 

25.7143 

12.5000 

0.51500 

GRID 

88 

25.7143 

12.5000 

0.50000 

GRID 

89 

25.7143 

10.0000 

0.51500 

GRID 

90 

25.7143 

10.0000 

0.50000 

GRID 

91 

25.7143 

7.50000 

0.51500 

GRID 

92 

25.7143 

7.50000 

0.50000 

GRID 

93 

25.7143 

5.00000 

0.51500 

GRID 

94 

25.7143 

5.00000 

0.50000 

GRID 

95 

25.7143 

2.50000 

0.51500 

GRID 

96 

25.7143 

2.50000 

0.50000 

GRID 

97 

25.7143 

0 . 

0.51500 

GRID 

98 

25.7143 

0 . 

0.50000 

GRID 

99 

22.5000 

15.0000 

0.51500 

G81D 

100 

22.5000 

15.0000 

0.50000 

GRID 

101 

22.5000 

12.5000 

0.51500 

GRID 

102 

22.5000 

12.5000 

0.50000 

GRID 

103 

2Z.5000 

10.0000 

0.51500 

GRID 

104 

22.5000 

10.0000 

0.50000 

GRID 

105 

22.5000 

7.50000 

0.51500 

GRID 

106 

22.5000 

7.50000 

0.50000 

GRID 

107 

22.5000 

5.00000 

0.51500 

GRID 

108 

22.5000 

5.00000 

0.50000 
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16.0714 ; 
16.0714 ; 
16.0714 l 
16.0714 I 
12.6571 . 
12.8571 , 
12.8571 ! 
12.8571 I 
12.8571 1 
12.8571 ] 
12.8571 1 
12.8571 ( 
12.8571 5 
12.8571 5 
12.8571 2 
12.8571 2 
12.8571 0 
12.8571 0 
9.64286 I 
9.64286 1 
6.42657 1 
6.42857 1 
3.21429 1 


»0 2.50000 0 
0 2.50000 0 
0 0. 0 
0 0 . 0 
7 15.0000 0 
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7 12.5000 0, 
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7 7.50000 0. 
7 7.50000 0, 
7 5.00000 0. 
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7 0. 0. 

7 0. 0. 
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I 10.0000 0. 
I 7.50000 0. 
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5.00000 0.5 
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0. 0.5 
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15.0000 0.5 
15.0000 0.5 
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GRID 

160 

3.21429 15.0000 0.50000 

GRID 

161 

0. 15.0000 0.51500 

GRID 

162 

0. 15.0000 0.50000 

GRID 

163 

9.64266 12.5000 0.51500 

GRID 

164 

9.64286 12.5000 0.50000 

GRID 

165 

9.64206 10.0000 0.51500 

GRID 

166 

9.64286 10.0000 0.50000 

GRID 

167 

9.64206 7.50000 0.51500 

GRID 

166 

9.64286 7.50000 0.50000 

GRID 

169 

9.64286 5.00000 0.51500 

GRID 

no 

9.64286 5.00000 0.50000 

GRID 

171 

9.64286 2.50000 0.51500 

GRID 

172 

9.64286 2.50000 0.50000 

GRID 

173 

9.64266 0. 0.51500 

GRID 

174 

9.64286 0. 0.50000 

GRID 

175 

6.42657 12.5000 0.51500 

GRID 

176 

6.42857 12.5000 0.50000 

GRID 

177 

3.21429 12.5000 0.51500 

GRID 

178 

3.21429 12.5000 0.50000 

GRID 

179 

0. 12.5000 0.51500 

GRID 

160 

0. 12.5000 0.50000 

GRID 

181 

6.42857 10.0000 0.51500 

GRID 

182 

6.42857 10.0000 0.50000 

GRID 

183 

6.42857 7.50000 0.51500 

GRID 

184 

6.42857 7.50000 0.50000 

GRID 

185 

6.42857 5.00000 0.51500 

GRID 

186 

6.42857 2.50000 0.51500 

GRID 

187 

6.42857 5.00000 0.50000 

GRID 

188 

6.42857 2.50000 0.50000 

GRID 

169 

6.42657 0. 0.51500 

GRID 

190 

6.42857 0. 0.50000 

GRID 

191 

3.21429 10.0000 0.51500 

GRID 

192 

3.21429 10.0000 0.50000 

GRID 

193 

0. 10.0000 0.51500 

GRID 

194 

0. 10.0000 0.50000 

GRID 

195 

3.21429 0. 0.51500 

GRID 

196 

3.2M29 0. 0.50000 

GRID 

J97 

0. 0. 0.51500 

GRID 

198 

0. 0. 0.50000 

GRID 

199 

3.214P9 7.50000 0.51500 

GRID 

200 

3.21429 5.00000 0.51500 

GRID 

201 

3.21429 2.50000 0.51500 

GRID 

202 

3.21429 7.50000 0.50000 

GRID 

203 

3.21429 5.00000 0.50000 

GRID 

204 

3.21429 2.50000 0.50000 

GRID 

205 

0. 7.50000 0.51500 

GRID 

206 

0. 5.00000 0.51500 

GRID 

207 

0. 2.50000 0.51500 

GRID 

208 

0. 7.50000 0.50000 

GRID 

209 

0. 2.50000 0.50000 

GRID 

210 

0. 5.00000 0.50000 


98 



CQiim 

1 

1 198 

C0UAD4 

2 

1 196 

CQUAD1 

3 

1 190 

C0UAD1 

4 

1 174 

CQUABI 

5 

1 154 

CO 1 -'API 

6 

1 140 

CQUAP1 

7 

1 126 

CQUAP4 

8 

1 112 

CQPAP4 

9 

1 90 

CQUAD4 

10 

1 84 

C0UAD1 

11 

1 70 

CQUAD4 

12 

1 56 

couAr* 

13 

1 42 

COUADA 

14 

1 28 

COUAD4 

15 

1 209 

C0UAD4 

16 

1 204 

CQUAD4 

17 

1 186 

C0UAD4 

18 

1 172 

C0UAD4 

19 

1 152 

COHAP4 

20 

1 138 

CODA D 4 

21 

1 124 

C0UAP4 

22 

1 110 

CQUAD4 

23 

1 96 

COUAD4 

24 

1 82 

CUUAD4 

25 

1 60 

COUAD4 

26 

1 54 

COUAD4 

27 

1 40 

C Q u A P 4 

28 

1 26 

CQUAD4 

29 

1 210 

CQUAD4 

30 

1 203 

CQUAP4 

31 

1 187 

C0UAD4 

32 

1 170 

CO'JAD4 

33 

1 150 

C0UAD4 

34 

1 136 

C0UAD4 

35 

1 122 

CQUAD4 

36 

1 108 

COUAD4 

37 

1 94 

CQUAD4 

38 

1 80 

CQIIAP4 

39 

1 66 

CQUAD4 

40 

1 52 

CQ1IAD4 

41 

1 38 

CQUAP4 

42 

1 24 

CQUAD4 

43 

1 200 

C 0 U A D 4 

44 

1 202 

COUAD4 

45 

1 184 

CODA D 4 

46 

1 168 

CQUAP4 

47 

1 148 

CQUAD4 

48 

1 134 

CQUAP4 

49 

1 120 

CQUAP1 

50 

1 106 

CQPAD4 

51 

1 92 


196 

204 

209 

-.25 

190 

188 

204 

-.25 

174 

172 

108 

-.25 

154 

152 

172 

-.25 

140 

130 

152 

-.25 

126 

124 

138 

-.25 

112 

110 

124 

-.25 

98 

96 

no 

-.25 

04 

62 

96 

-.25 

70 

68 

82 

-.25 

56 

54 

60 

-.25 

42 

40 

54 

-.25 

28 

25 

4C 

-.25 

14 

12 

26 

-.25 

204 

203 

210 

-.25 

188 

187 

203 

-.25 

172 

170 

187 

-.25 

152 

150 

170 

-.25 

130 

136 

150 

-.25 

124 

122 

136 

-.25 

no 

108 

122 

-.25 

96 

94 

108 

-.25 

82 

60 

94 

-.25 

68 

66 

80 

-.25 

54 

52 

66 

-.25 

40 

38 

52 

-.25 

26 

24 

38 

-.25 

12 

10 

24 

-.25 

203 

202 

208 

-.25 

187 

184 

202 

-.25 

170 

168 

184 

-.25 

150 

148 

168 

-.25 

136 

13* 

148 

-.25 

122 

120 

134 

-.25 

108 

106 

120 

-.25 

94 

92 

106 

-.25 

00 

78 

92 

-.25 

66 

64 

70 

-.25 

52 

50 

64 

-.25 

38 

36 

50 

-.25 

24 

22 

36 

-.25 

10 

8 

22 

-.25 

202 

192 

194 

-.25 

184 

182 

192 

-.25 

168 

166 

102 

-.25 

140 

146 

166 

-.25 

134 

132 

146 

-.25 

120 

118 

132 

.25 

106 

104 

118 

-.25 

92 

90 

104 

-.25 

70 

76 

90 

-.25 


99 




1 


CQUAD4 

CQIIAP4 

CQUAP4 

CCIIAD4 

CO PAD 4 

COUAD4 

COOAD4 

COUAP4 

CQUAD4 

CQUAD4 

CQHAD4 

CQUAD4 

CQUAD4 

CQUAD4 

CQUAD4 

CQUAD4 

COUAP4 

CQUAD4 

CQUAD4 

CQUAD4 

COUAD4 

CQUAD4 

CQUAD4 

CQUAD4 

CQUAD4 

COUAD4 

CQHAD4 

COUAD4 

C0UAD4 

CQUAD4 

CQUAD4 

CQUA04 

CQDAD4 

CQUAD4 

CQUAD4 

CQUAD4 

COUAD4 

CQUAD4 

CQUAD4 

CQUAD4 

COUAD4 

COUAD4 

CQUAD4 

COHAD4 

CQUAD4 

CQ0AD4 

COUAD4 

COWAD4 

COUAD4 

COHAD4 

COUAD4 


62 
63 

54 

55 

56 i 

57 1 

58 1 

59 J 

60 1 

61 1 

62 1 

63 l 

64 1 

65 1 

66 1 

67 1 

68 1 

69 1 

70 1 

71 1 

72 1 

73 1 

74 1 

75 1 

76 1 

77 1 

78 . 1 

79 1 

80 1 

81 1 

82 1 

83 1 

84 1 

85 85 

86 85 

87 85 

88 85 

89 85 

90 65 

91 85 

92 85 

93 85 

94 85 

95 85 

96 85 

97 85 

90 85 

99 85 

100 85 

101 85 

102 85 


64 

50 

50 

36 

36 

22 

22 

8 

194 

192 

192 

182 

182 

166 

166 

146 

116 

132 

132 

118 

118 

IS4 

104 

90 

90 

76 

76 

62 

62 

40 

48 

34 

34 

20 

20 

6 

180 

178 

178 

176 

176 

164 

164 

144 

144 

130 

130 

116 

116 

102 

102 

88 

88 

74 

74 

60 

60 

46 

46 

32 

32 

18 

18 

4 

197 

195 

195 

189 

189 

173 

173 

153 

153 

139 

139 

125 

125 

111 

HI 

97 

97 

83 

83 

69 

69 

55 

55 

41 

41 

27 

27 

13 

207 

201 

201 

186 

186 

17! 

171 

151 


48 62 

34 48 

20 34 

6 20 

178 100 

176 178 

164 176 

144 164 

130 144 

116 130 

102 116 

88 102 

74 88 

60 74 

46 60 

32 46 

18 32 

4 16 

160 162 

158 160 

156 150 

142 156 

128 142 

114 128 

100 114 

86 100 

72 86 

56 72 

44 58 

30 44 

16 30 

2 16 

201 207 

186 201 

171 186 

151 171 

137 151 

123 137 

109 123 

95 109 

81 95 

67 81 

53 67 

39 53 

25 39 

11 25 

200 206 

185 200 

169 185 

149 169 


-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

.25 

-.25 

-.25 

-.25 

-.25 

-.25 

-.25 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 


1 00 





f :OllAD< 

CQUAD4 

CQUAD4 


103 

10 < 

105 













































C0UAD4 

154 

CQUAD1 

155 

CQUAD4 

156 

CQUAD4 

157 

C0UAD4 

158 

CQUAD4 

159 

CQUAD4 

160 

CQUAD4 

161 

CQ11AD4 

162 

CQ0AD4 

163 

COOAD4 

164 

CQUAD4 

165 

CQUAD4 

166 

C0UAD4 

167 

CQUAD4 

168 

CHEIA 

169 

<E 169 

201 

CHEIA 

170 

«E 170 

186 

CHEIA 

171 

*E 171 

171 

CHEIA 

172 

*E 172 

151 

CHEIA 

173 

«E 173 

137 

CHEIA 

174 

• E 174 

123 

CHEIA 

175 

<E 175 

109 

CHEIA 

176 

*E 176 

75 

CHEIA 

177 

«E 177 

81 

CHEIA 

178 

»E 178 

67 

CHEIA 

179 

*E 179 

53 

CHEIA 

180 

♦ E 180 

39 

CHEIA 

181 

•E 181 

25 

CHEIA 

102 

•E 102 

11 

CHEIA 

103 

*E 103 

200 

CHEIA 

184 

*E 164 

185 

CHEIA 

185 

CD 

cn 

169 

CHEIA 

186 

•E 186 

149 


85 

05 

85 

05 

85 

05 

85 

85 

85 

05 

85 

05 

85 

85 

85 

2 

207 

2 

201 

2 

186 

2 

171 

2 

151 

2 

137 

■2 

123 

2 

109 

2 

95 

2 

81 

2 

67 

2 

53 

2 

39 

2 

25 

2 

206 

2 

200 

2 

185 

2 

169 


19 

179 

177 

175 

163 

H3 

129 

115 

101 

07 

73 

59 

45 

31 

17 


196 

190 

174 

154 

140 

126 

112 

98 

84 

70 

56 

42 

28 

209 

204 

188 

172 


5 

3 

17 

177 

159 

161 

175 

157 

159 

163 

155 

157 

143 

141 

155 

129 

127 

111 

115 

113 

127 

101 

99 

113 

07 

85 

99 

73 

71 

85 

59 

57 

71 

45 

43 

57 

31 

29 

43 

17 

15 

29 

3 

1 

15 

196 

204 

209 

190 

180 

204 

174 

172 

188 

154 

152 

172 

140 

138 

152 

126 

124 

138 

112 

no 

124 

96 

96 

110 

84 

82 

96 

70 

68 

82 

56 

54 

60 

42 

40 

54 

20 

26 

40 

14 

12 

26 

204 

203 

210 

188 

187 

203 

172 

170 

187 

152 

150 

170 


.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 

.125 


197 

195 

E 

169 

195 

189 

E 

170 

189 

173 

E 

171 

173 

153 

E 

172 

153 

139 

E 

173 

139 

125 

E 

174 

125 

111 

E 

175 

111 

97 

E 

176 

97 

83 

E 

177 

03 

69 

E 

178 

69 

55 

E 

179 

55 

41 

E 

180 

41 

27 

E 

161 

27 

13 

E 

182 

207 

201 

E 

183 

201 

106 

E 

184 

186 

171 

E 

185 

171 

151 

E 

186 


102 



CHEIA 

ie7 

2 

152 

138 

t£ 187 

135 

149 



CHEKA 

188 

2 

138 

124 

• £ 188 

121 

135 



CHEKA 

109 

2 

124 

no 

• E 189 

107 

121 



CHEKA 

190 

2 

110 

96 

•E 190 

93 

107 



CHEKA 

191 

2 

96 

82 

«E 191 

79 

93 



CHEIA 

192 

2 

82 

68 

♦ F 192 

65 

79 



CHEIA 

193 

2 

68 

54 

♦ E 193 

51 

65 



CHEIA 

194 

2 

54 

40 

• E 194 

37 

51 



CHEIA 

195 

2 

40 

26 

♦ E 195 

23 

37 



CHEIA 

196 

2 

26 

12 

«E 196 

9 

23 



CHEIA 

197 

2 

210 

203 

♦ E 197 

199 

205 



CHEIA 

198 

2 

203 

187 

«E 198 

183 

199 



CHEIA 

199 

2 

187 

170 

♦ E 199 

167 

183 



CHEIA 

200 

2 

no 

ISO 

♦ E 200 

147 

167 



CHEIA 

201 

2 

150 

136 

»E 201 

133 

147 



CHEIA 

202 

2 

136 

122 

*E 202 

119 

133 



CHEIA 

203 

2 

122 

108 

♦ E 203 

105 

119 



CHEIA 

204 

2 

108 

94 

♦ E 204 

91 

105 



CHEIA 

205 

2 

94 

80 

♦ E 205 

77 

91 



CHEIA 

206 

2 

80 

66 

•E 206 

63 

77 



CHEIA 

207 

2 

66 

52 

•E 207 

49 

63 



CHEIA 

208 

2 

52 

38 

«E 208 

35 

49 



CHEIA 

209 

2 

38 

24 

♦ E 209 

21 

35 



CHEIA 

210 

2 

24 

10 

♦E 210 

7 

21 



CHEIA 

211 

2 

208 

202 

'£ 211 

19* 

193 




136 

150 

151 

137 

E 

187 

122 

136 

137 

123 

E 

108 

108 

122 

123 

109 

E 

189 

94 

108 

109 

95 

E 

190 

80 

94 

95 

81 

E 

191 

66 

80 

81 

67 

E 

192 

52 

66 

67 

53 

E 

193 

38 

52 

53 

39 

E 

194 

24 

38 

39 

25 

E 

195 

10 

24 

25 

11 

E 

196 

202 

208 

206 

200 

E 

197 

164 

202 

200 

185 

E 

198 

168 

184 

185 

169 

E 

199 

148 

168 

169 

149 

E 

200 

134 

146 

149 

135 

E 

201 

120 

134 

135 

121 

E 

202 

106 

120 

121 

107 

E 

203 

92 

106 

107 

93 

E 

204 

78 

92 

93 

79 

E 

205 

64 

70 

79 

65 

E 

206 

50 

64 

65 

51 

E 

207 

36 

50 

51 

37 

E 

208 

22 

36 

37 

23 

E 

209 

8 

22 

23 

9 

E 

210 

192 

194 

205 

199 

E 

2)1 


103 



CHEKA 

212 

2 

202 

184 

<E 212 

181 

191 



CHEKA 

213 

2 

184 

168 

*E 213 

165 

161 



CHEKA 

214 

2 

168 

148 

iE 2H 

145 

165 



CHEKA 

215 

2 

140 

134 

*E 215 

131 

145 



CHEKA 

216 

2 

134 

120 

• E 216 

117 

131 



CHEKA 

217 

2 

120 

106 

'£ 217 

103 

117 



CHEKA 

218 

2 

106 

92 

«£ 212 

89 

103 



CHEKA 

219 

2 

92 

78 

«E 219 

75 

89 



CHEKA 

220 

2 

78 

64 

♦ E 220 

61 

75 



CHEKA 

221 

2 

64 

50 

«E 221 

47 

61 



CHEKA 

222 

2 

50 

36 

»E 222 

33 

47 



CHEKA 

223 

2 

36 

22 

«E 223 

19 

33 



CHEKA 

224 

2 

22 

8 

♦ E 224 

5 

19 



CHEKA 

225 

2 

194 

192 

• E 225 

177 

179 



CHEKA 

226 

. 2 

192 

182 

*E 226 

175 

177 



CHEKA 

227 

2 

182 

166 

«E 227 

163 

175 



CHEIA 

228 

2 

166 

146 

»B 228 

143 

163 



CHEIA 

229 

2 

146 

132 

*E 229 

129 

143 



CHEIA 

230 

2 

132 

lie 

«E 230 

115 

129 



CHEIA 

231 

2 

118 

104 

*E 231 

101 

115 



CHEKA 

232 

2 

104 

90 

•E 232 

87 

101 



CHEIA 

233 

2 

90 

76 

♦ E 233 

73 

87 



CHEIA 

234 

2 

76 

62 

«E 234 

59 

73 



CHEIA 

235 

2 

62 

48 

♦ E 235 

45 

59 



CHEIA 

236 

2 

48 

34 

<E 236 

31 

45 



CHEKA 

237 

2 

34 

20 

• E 237 

17 

31 




182 

192 

199 

183 

E 

212 

166 

182 

183 

167 

E 

213 

146 

166 

167 

147 

E 

214 

132 

146 

147 

133 

E 

215 

118 

132 

133 

119 

E 

216 

104 

118 

119 

105 

E 

217 

90 

104 

105 

91 

E 

218 

76 

90 

91 

77 

E 

219 

62 

76 

77 

63 

E 

220 

48 

62 

63 

49 

E 

221 

34 

48 

49 

35 

E 

222 

20 

34 

35 

21 

E 

223 

6 

20 

21 

7 

E 

224 

178 

160 

193 

191 

E 

225 

176 

176 

191 

181 

E 

226 

164 

176 

161 

165 

E 

227 

144 

164 

165 

145 

E 

228 

130 

144 

145 

131 

E 

229 

116 

130 

131 

117 

E 

230 

102 

116 

117 

103 

E 

231 

88 

102 

103 

89 

E 

232 

74 

88 

89 

75 

E 

233 

60 

74 

75 

61 

E 

234 

46 

60 

61 

47 

E 

235 

32 

46 

47 

33 

E 

236 

18 

32 

33 

19 

E 

237 


104 





CBEXA 

230 

2 

20 

6 

4 

18 

19 

5 E 

238 

<E 238 

3 

17 








CREXA 

239 

2 

180 

178 

160 

162 

179 

177 E 

239 

»E 239 

159 

161 








CHEXA 

240 

2 

178 

176 

158 

160 

177 

175 E 

240 

• E 240 

157 

159 








CflEIA 

241 

2 

176 

164 

156 

158 

175 

163 E 

241 

•E 241 

155 

157 








CHEXA 

242 

2 

164 

144 

142 

156 

163 

143 E 

242 

<E 242 

141 

155 








CHEXA 

243 

2 

144 

130 

120 

142 

143 

129 E 

243 

• £ 243 

127 

141 








CHEXA 

2*4 

2 

130 

1 16 

114 

128 

129 

115 E 

244 

*E 244 

113 

127 








CHEXA 

245 

2 

116 

102 

100 

114 

115 

101 E 

245 

♦E 245 

99 

113 








CHEXA 

246 

2 

102 

88 

86 

100 

101 

87 E 

246 

<E 246 

85 

99 








CHEXA 

247 

2 

80 

74 

72 

86 

87 

73 E 

247 

*E 247 

71 

85 








CHEXA 

248 

2 

74 

60 

58 

72 

73 

59 E 

248 

• E 240 

57 

71 








CHEXA 

249 

2 

60 

46 

44 

58 

59 

45 E 

249 

• E 249 

43 

57 








CHEXA 

250 

2 

46 

32 

30 

44 

45 

31 E 

250 

»E 250 

29 

43 








CHEXA 

251 

2 

32 

18 

16 

30 

31 

17 E 

251 

*E 251 

15 

29 








CHEIA 

252 

• 2 

18 

4 

2 

16 

17 

3 E 

252 

*E 252 

1 

15 








PSHELL 

1 

1 0.50000 

1 

1.12200 

1 




PSHELL 

85 

1 0.25000 

1 

1.12200 

1 




psolld,2,2 










«atl,l,1.0e7, 

,0.33,Z.58e-4 








■atl,2,,251.4 

1,0.49,9.067e-5 








elgr,l,»giv,C 

1.0,5000.0,,30,, 

,’elgr 








•elgr.iass 
suport,106,12345 
3 pci,1,6,1,thru,210 
parai.autospc.yes 
paralysing,1 
rloadl,10,2001,,,3001 
da rea,2001,79,3,1.0 

tabledl ,3001.,«tal 

♦ tal, 0.0,1.0,5000.0,1.0,endt 
(reql,10,5.0,1.0,1000 
tabd*pl,101,,,,,,,,'d»p 

*d»p,40.0,0.214,63.0,0.216,95.0,0.218,160.0,0.222,«d»pl 
«d»pl,200.0,0.224,285.0,0.225,320.0,0.224,440.0,0.216,«d»p2 
«d*p2,460.0,0.215,480.0,0.212,520.0,0.208,580.0,0.197,<dip3 
*uip3,630.0,0.109,665.0,0.163,715.0,0.171,825.0,0.142,«d«p4 
•d*p4,870.0,0.130,885.0,0.125,1000.0,0.085,endt 


105 










LIST OF REFERENCES 


1. Johnson, C.D. and Kienholz, D.A., “Finite Element Prediction of 
Damping in Structures With Constrained Viscoelastic Layers,” AIAA 
Journal, Vol. 20, No. 9, September 1982. 

2. Maurer, G.J., Vibration R esponse of Constrained Viscoelasticallv 
Damped Plates: Analysis and Experiments . M.S. Thesis, U.S. Navai 
Postgraduate School, Monterey, California, December 1987. 

3. Nashif, A.D., Jones, D.I.G., and Henderson, J.P., Vibration Damping . 
Wiley-Interscience Publications, John Wiley and Sons, 1985. 

4. Jones, D.I.G., and Henderson, J.P., “Fundamentals of Damping 
Materials,” Vibration Damping Short Course Notes, Section 2.3, 
University of Dayton, Dayton, Ohio, June 1987. 

5. Rogers, L., “Single Constrained Layer Damping Treatment Analysis,” 
Vibration Damping Short Course Notes, Section 6.3, University of 
Dayton, Dayton, Ohio, June 1987. 

6. Henderson, J.P., and Jones, D.I.G., “Fundamentals of Layered 
Damping Treatments,” Vibration Damping Short Course Notes, 

Section 6.1, University of Dayton, Dayton, Ohio, June 1987. 

7. MacNeal-Schwendler Corporation, MSC/NASTRAN Users Manual . 
Vol. 1, MacNeal-Schwendler Corporation, November, 1985. 

8. Johnson, C.D., and Kienholz, D.A., “Finite Element Design of 
Viscoelastically Damped Structures,” CSA Engineering, Inc., Palo 
Alto, California, 1984. 

9. Nashif, A.D., “Layered Damping Treatment Design, Multiple Layer,” 
Vibration Damping Short Course Notes, Section 6.4, University of 
Dayton, Dayton, Ohio, June 1987. 

10. Ross, D., Ungar, E.E., and Kerwin, E.M., “Damping of Plate Flexural 
Vibrations by Means of Viscoelastic Laminae,” Structural Damping . 
ASME, 1959. 

11. Drake, M.L., “Design Techniques - Fourth Order Beam Theory,” 
Vibration Damping Short Course Notes, Section 7.1, University of 
Dayton, Dayton, Ohio, 1987. 


106 






12. University of Dayton Research Institute, “ISD-112 Material Data 
Sheet,” Vibration Damping Short Course Notes, University of Dayton, 
Dayton, Ohio, June 1987. 

13. Gere, J.M., and Timoshenko, S.P., Mechanics of Materials . 2nd 
Edition, PWS Publishers, Boston, Massachusetts, 1984. 

14. Richardson, N., and Potter, R., “Identification of the Modal Properties 
of an Elastic Structure From Measured Transfer Function Data,” 
presented at the 20th International Instrumentation Symposium, May 
1984, Albuquerque, New Mexico, Instrument Society of America 
reprint, 1974. 

15. PDA Engineering, PAT/MSC Interface Guide . Release 2.1, PDA 
Engineering, September 1988. 

16. Johnson, C.D., and Kienholz, D.A., “Finite Element Prediction of 
Damping in Structures With Constrained Viscoelastic Layers,” Report 
No. 82.053, Anamet Laboratories, May 1982. 

17. MacNeal-Schwendler Corporation, MSC/NASTRAN Handbook for 
Dynamic Analysis . Section 3.2.3, MacNeal-Schwendler Corporation, 
1983. 

18. PDA Engineering, PATRAN PLUS Users Manual . Vol. 2, Ch. 21, PDA 
Engineering, 1987. 


107 








INITIAL DISTRIBUTION LIST 


•» 

1. 

Defense Technical Information Center 

Cameron Station 

Alexandria, Virginia 22304-6145 

2 

i 

2. 

Library, Code 0142 

Naval Postgraduate School 

Monterey, California 93943-5002 

2 


3. 

Dean of Science and Engineering, Code 06 

Naval Postgraduate School 

Monterey, California 93943-5004 

2 


4. 

Research Administrations Office, Code 012 

Naval Postgraduate School 

Monterey, California 93943-5004 

1 


5. 

Department Chairman, Code 69 

Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, California 93943-5004 

1 

• 

6. 

Curricular Officer, Code 34 

Naval Engineering Curriculum 

Naval Postgraduate School 

Monterey, California 93943—4004 

1 

f 

7. 

Professor Y.S. Shin, Code 69Sg 

Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, California 93943-5004 

3 


8. 

Dr. K.S. Kim, Code 69Ki 

Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, California 93943-5004 

1 


9. 

Dr. Arthur Kilcullen, Code 1962 

David W. Taylor Naval Ship R&D Center 

1 



Bethesda, Maryland 20084 


* 


108 











10. Dr. Lawrence Maga, Code 196 2 

David W. Taylor Naval Ship R&D Center 

Bethesda, Maryland 20084 

11. Mr. Doug Noll, Code 1944 1 

David W. Taylor Naval Ship R&D Center 

Bethesda, Maryland 20084 

12. LT Michael J. Bateman, USN 1 

5306 SE 64th 

Portland, Oregon 97206 


109 



